Volume 19, No. 2 


ABTOMATUKA 
February 1958 
TEAEMEXAHHKUKA 
OF 
MAY 4 - 1969 
SOVIET INSTRUMENTATION AND 


CONTROL TRANSLATION SERIES 


Automation 


Control 


(The Soviet Journal Avtomatika i Telemekhanika in English Translation) 


@ This translation of a Soviet journal on automatic control is 
published as a service to American science and industry. It is 
sponsored by the Instrument Society of America under a grant 
in aid from the National Science Foundation, continuing a pro- 
gram initiated by the Massachusetts Institute of Technology. 


' 
2 
i 
' 
4 
- 
q 
\ 
q 


SOVIET INSTRUMENTATION AND 
CONTROL TRANSLATION SERIES 


Instrument Society 
of America 
Executive Board 


Henry C. Frost 
President 

Robert J. Jeffries 

Past President 

John Johnston, Jr. 
President-Elect-Secretary 
Glen G. Gallagher 
Dept. Vice President 
Thomas C. Wherry 

Dept. Vice President 
Philip A. Sprague 
Dept. Vice President 
Ralph H. Tripp 

Dept. Vice President 
Howard W. Hudson 
Treasurer 

Willard A. Kates 
Executive Assistant—Districts 
Benjamin W. Thomas 
Executive Assistant—Conferences 
Carl W. Gram, Jr. 
Dist. I Vice President 
Charles A. Kohr 

Dist. II Vice President 
J. Thomas Elder 

Dist. III Vice President 
George L. Kellner 
Dist. IV Vice President 
Gordon D. Carnegie 
Dist. V Vice President 
Glenn F. Brockett 
Dist. VI Vice President 
John F. Draffen 

Dist. VII Vice President 


John A. See 

Dist. VIII Vice President 
Adelbert nter 
Dist. IX Vice gag 
Joseph R. 

Dist. X Vice roident 
Headquarters Office 


William H. Kushnick 
Executive Director 


Charles W. Covey 

Editor, ISA Journal 

George A. Hall, Jr. 

Assistant Editor, ISA Journal 
Herbert S. Kindler 

Director, Tech. & Educ. Services 
Ralph M. Stotsenburg 
Director, Promotional Services 


William F. Minnick, Jr. 
Promotion Manager 


ISA Publications Committee 


Nathan Cohn, Chairman 
Jere E. Brophy Richard W. Jones John E. Read 
Enoch J. Durbin George A. Larsen Joshua Stern 


George R. Feeley Thomas G. MacAnespie Frank S. Swaney 
Richard A. Terry 


Translations Advisory Board 
of the Publications Committee 


Jere E. Brophy, Chairman 
T. J. Higgins S. G. Eskin G. Werbizky 


@ This translation of the Soviet Journal Avtomatika i 
Telemekhanika is published and distributed at nominal 
subscription rates under a grant in aid to the Instrument 
Society of America from the National Science Foundation. 
This translated journal, and others in the Series (see back 
cover), will enable American scientists and engineers to 
be informed of work in the fields of instrumentation, 
measurement techniques and automatic control reported 
in the Soviet Union. 


The original Russian articles are translated by competent 
technical personnel. The translations are on a cover-to- 
cover basis, permitting readers to appraise for themselves 
the scope, status and importance of the Soviet work. 


Publication of Avtomatika i Telemekhanika in English 
translation started under the present auspices in April 
1958 with Russian Vol. 18, No. 1 of January 1957. Trans- 
lation of Vol. 18 has now been completed. The twelve issues 
of Vol. 19 will be published in English translation by 
mid-1959. 


All views expressed in the translated material are intended 
to be those of the original authors, and not those of the 
translators, nor the Instrument Society of America. 


Readers are invited to submit communications on the qual- 
ity of the translations and the content of the articles to 
ISA headquarters. Pertinent correspondence will be pub- 
lished in the “Letters” section of the ISA Journal. Space 
will also be made available in the ISA Journal for such 
replies as may be received from Russian authors to com- 
ments or questions by American readers. 


Subscription Prices: 
Per year (12 issues), starting with Vol. 19, No. 1 


General: United States and Canada. . . . . $30.00 
Libraries of non-profit institutions: 
United States and Canada. . . . . $15.00 
Single issues to everyone, 


See back cover for combined subscription to entire Series. 


Subscriptions and requests for information on back issues 
should be addressed to the: 

Instrument Society of America 

313 Sixth Avenue, Pittsburgh 22, Penna. 


Translated and printed by Consultants Bureau, Inc. 


Volume XIX No.2 — February 1958 


English Translation Printed January 1959 


Automation and Remote Control 


The Soviet Journal Avtomatika i Telemekhanika 


in English Translation 


Avtomatika i Telemekhanika is a Publication of the Academy of Sciences of the USSR 


EDITORIAL BOARD 
as Listed in the Original Soviet Journal 


Corr. Mem. Acad. Sci. USSR V. A. Trapeznikov, Editor in Chief 
Dr. Phys. Math. Sci. A. M. Letov, Assoc. Editor 

Academician M. P. Kostenko 

Academician V. S. Kulebakin 

Corr. Mem. Acad. Sci. USSR B. N. Petrov 


Dr. 
Dr. 
Dr. 
Dr. 
Dr. 


Dr 


Tech. Sci. M. A. Aizerman 
Tech. Sci. V. A. Il’in 

Tech. Sci. V. V. Solodovnikov 
Tech. Sci. B. S. Sotskov 
Tech. Sci. Ia. Z. Tsypkin 
Tech. Sci. N. N. Shumilovskii 


Cand. Tech. Sci. V. V. Karibskii 

Cand. Tech. Sci. G. M. Ulanov, Corresp. Secretary 
Eng. S. P. Krasivskii 

Eng. L. A. Charikhov 


See following page for Table of Contents. 
Copyright by Instrument Society of America 1959 


5 
i 
4 
| 
q 
f 


AUTOMATION AND REMOTE CONTROL 


VOLUME 19, NUMBER 2 FEBRUARY 1958 


CONTENTS 


On The Stability of Nonlinear Control Systems Described by Fifth and Sixth Order 
Differential Equations, E, N. Rozenvasser 91 


On The Theory of One-Relay Systems. V. 5S, Boiarinov and N, N, Leonov ......+6++++ 103 


Determination of Parameters for Corrective Devices in Linear Servosystems Using 
Given Generalized Parameters. M, M, Kreimerman ...... 


124 


An Analytical Forrnulation of the Synthesis Problem of Corrective Devices in 
Linear Servosystems, V. G, Segalin eee 137 


An AC Bridge With Self-Balancing In Two Parameters, V. lu, Kneller ........++++- 156 


Frequency Methods of Remote Control For Distributed Objects, V. A, Il'in and 
K, P, Kurdiukov 167 


Letters to the Editor 


A Nonlinear Semiconductor of Resistance Sensitive to Magnetic Fields, G, E, Pikus 
and O, V. Sorokin 177 


Chronicle 


The International Federation of Automatic Control (IFAC), A, M, Letov and 

B, N, Naumov 179 
Conference on Automatic Control and Computer Technique. I, M, Rusevich  ........ 182 
Errata 


186 


162 


174 


187 


189 
191 
195 


| 
RUSS. 
PAGE PAGE 
114 
fr 
4 the 
cit 
of | 
| 
} 


ON THE STABILITY OF NONLINEAR CONTROL SYSTEMS DESCRIBED BY 
FIFTH AND SIXTH ORDER DIFFERENTIAL EQUATIONS 


E, N, Rozenvasser 


(Leningrad) 


On the basis of a theorem due to Lure [1], sufficiency conditions “in the large” 
are obtained for some control systems described by fifth- and sixth-order differential 
equations, 


Let the control systems under consideration be described by the Equations {1,1,19} * 


= + haf (k= 1, 


n (1) 
c= (nm = 5, 6). 
s=1 
We assume that the characteristic equation of the open System (1) 
D (i) = | bea — na | = 0 (2) 


lfork= a 


(5 key is the Kronecker symbol, & 


from zero and with negative real parts,** 
The characteristic equation of the linearized System (1), with f(g) =co, may be put in the form 


4 (u) = — cM (p), (3) 


where the degree of the polynomial M(y) is always less than the degree of the polynomial D(y), In the case 
n=6, it will be assumed that 


) has roots A=Ap =1,2,..., 0) (n=5, 6) different 


(4) 
This relation will hold identically if the degree of M(u) does not exceed three, For simplicity we assume 
that this is indeed the case, Thus, we may assume: 
for n=5 


D (p) = — + — + — + key, (5) 
M (p) = — + — + — (6) 


* Formulas and notation taken from Lur'e's book (1) will retain the numeration therein give... The corresponding j 
Citations are given within braces, the Roman numeral giving the chapter, the Arabic numerals giving the number 
of the formula within the chapter, ' 


** The case will also be considered where n=6, \g=0,r=— Ljshs =0. 
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for n=6 


= — myp + — mau? + — + Me, 
M (yp) = my? — + — my. (8) 


According to a theorem of Lur'e [1], sufficiency conditions for stability “in the large” for systems (1) 
can be obtained as necessary and sufficient conditions for the existence of real solutions of a certain system of 
quadratic equations, This latter system is called the resolving system. In considering the case corresponding to 
{11.8, 16-21} , and taking (4) into account, we may put the resolving systems in the form: 


for n=5 


(to4)* + r, =0, 
+ — + + + T's = 0, 
— (te)? +P =0, 
+ — + + +T_s = 0; 
for n=6 
(t:9,)? + Rs =0, 
— (ty%q + + 2t + + + Ry =0, 
— +R, =9, 


+ 49_1)* — (tg9_3 + + + R_»=0, (10) 
— + + — (tg9_g + + + + + 
+ (t59_4 + + + + =0. 


In Systems (9) and (10) the tj are the unknowns, 


Lo, 1304, r; L129 + 1,9, 159, 
Dy = — + — + 1435 — 1595, 


Ry = Ry = 1,9, — + — 1435, 
Ry = R.; == — Mg + — (12) 


In what follows, we shall assume that [';, 3, Rs and R_y are all different from zero; otherwise we would 
be back to the case previously investigated[1]. The constants 0,(k=1, 2,..., 0+ 3), O.,(k = 1, 2, 3, 4) are 
computed from the coefficients of (5) and (7) in accordance with {II,8,12-13} , and the quantities o_, and 
Gg , obtained by methods given in [1], equal 


Sy = Ongy = — ki + Shik, — — + + — ks. (14) 


For the existence of real solutions of the resolving systems (9) and (10), it is obviously necessary that the 
following inequalities hold: 


for n=5 
r, > 0, r > 0, (15) 


and | 


for n =6 


Rs<0, R,>0. (16) 
In this case we have for System (9) 
and for System (10) 


If (17) and (18) hold, then systems (9) and (10) may be transformed to the single form: 


y—ez=B, 2#—2y=C, (19) 
with the following notation used: 
forn =5 
2V 2| to|V [total [total 4V | tote| (20) 
(2046, — 02) — T's ti (20_10_5 ,) 
Cm + 8A | tol tate | — 0G) — (05 + 20405 — + Ts cm) 
160% | tot, | : 
for n =6 
| 2V 2) | 4V| tits (22) 


Aw — + (20_,0_, —o1,) — R_, 
803 | ts | tits | 


C = 2tytyo* + 8B | ts | | tats | (26_,0_,—o” ,)— + R_, (23) 
1607, | tits| 


Taking the various combinations of signs in (19), we obtain four systems of quadratic equations, only two 
of which are essentially different, for example: 


2—2zy=C (24) 


a—z=A, y+2z=B, ay =C. (25) 


The other two systems are obtained from (24) and (25) by replacing z by ~2z, 


We will now show that, if for some A, B, and C, System (25) has real solutions, then System (24) also has 
such solutions, Considering System (25), one easily finds that, for the existence of real solutions for it, it is 
necessary and sufficient that the following relationship holds 


(2s) 
and that the equation 
$1 (z) = (2? — C)? — (A + 2) (B—z) = 0, (27) 
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obtained by eliminating x and y from System (25), has real roots. 
Indeed, if z is a real root of Equation (27), then 


(z, + A)(B—2) >0. (28) 


If now the direction of the inequality is the same as in the Condition (26), then one of the quantities x” 
and y* is positive (for example, x’ > 0). Then we obtain from (25), z1 + A > 0, and from (28), — y*=(B— 2) = 
= 0, i.e., the root z indeed provides a real solution to the system. If A+B =0, then by analogous reasoning 
one easily shows that there is a real solution only if A=B=C=0, Therefore, we shall henceforth assume that 
(26) is always satisfied. A result of the previous proof is that, given (26), any real root of Equation (27) gives a 
real solution to (25). 


If A and B are interchanged in (25), Equation (27) can be written as 
$2 (2) = (2? —C)* — (B z)(A—2) =0 (29) 


and will have roots with the opposite signs of those of Equation (27), 


It is obvious from this that the newly obtained system also will have real solutions, Consequently, if (25) 
has real solutions, then one of the Equations (27) or (29) will have nonnegative roots, Now, let us return to (24), 
It has real solutions if and only if the equation 


(2) = (2? —C)? — (2 + A) (2+ B) = 
has a real root 2, for which 
%+AD>O0 (31) 
and 
2+B>0. (32) 


If one of these inequalities is satisfied then, by virtue of (30), so is the other. We now prove that if (25) 
has real solutions, then Equation (30) has nonnegative roots, From (28), (29), and (30), we have 


? (2) = 9, (2) — 22(A + 2), (33) 
= — 22(B + 2). (34) 
From (26), there are two possible cases: 


A>0, B>0, (35) 
A>0, B<0. (36) 


The cases A > 0, B = 0, and A < 0, B> 0 are analogous to (35) and (36), If (35) holds and 2 is a non- 
negative root of Equation (27), then from (33) we have 


(29) = — 22 (A + 2) <0. 


Since yg a) > 0, Equation (30) in this case has a nonnegative root, zg = Zo. If (36) holds, then Equation 
(27) has only negative roots, and Equation (29) has only positive roots. if z's is one of the roots of Equation (27), 
then 


B+ 2420, (37) 


since, as was proved, if (26) holds then any real root of Equation (27) gives real solutions to (25), 
Substituting the value of z', in (34), we get 


A 
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that is, Equation (30) has a positive root, zg = z'». Whence, from (26), follows the existence of real solutions 
of (24), since this implies that the sum of A or B with the nonnegative root will be positive, for example, 
> 0 and, by virtue of (30), B+z, = 0. 


Thus, in order that the solution of the resolving systems, (9) and (10), be real, given (15) and (16), it is 
necessary and sufficient that (24) have real solutions. Thus, solving the problem of the stability of the trivial 
solution of (1) leads to the investigation of the quadratic Equations (24), From a detailed consideration of 
(24) we formulate our final result as a theorem, to be proved below, 


Let either (15) or (16) hold, Then, if 
C++4+AB>0 


+ ace — — 
— 27 | 2 _ +O — AP (44 BY|>0, (39) 


then for stability “in the large" of the unperturbed motion corresponding to (1), it is sufficient that at least one 
of the following conditions does not hold: 


A+B<0, AB>C>4; 


(40) 
if 
D<0 (41) 
or 
D=0 (42) 
and the following relationships do not hold 
C+3<0, C+i+4B=0, 44+B=0, (43) 
then it is sufficient for stability that the following inequality holds: 
A+B>0. (44) 


Thus, if (42) holds then (38) does not, The stability criteria obtained are the most general of all those 
that might be obtained for the system under consideration by the use of the Lur’e theorem, 


Proof, Consider (24), Eliminating x and y we obtain Equation (30), which can be put in the form: 
z*— (2C 4+ 1)2*— (A+ B)z4+C*— AB=0. (45) 


Using known rules [2], we find that if (38) and (39) hold, then Equation (45) has four real roots, If one 
of the Conditions (41) or (42) of the theorem to be proved holds, then Equation (30) has two real roots. To 
solve the problem posed, it is necessary to clarify whether or not (31) and (32) hold for at least one of the real 
roots of Equation (45) for each A and B, 


We make the following preliminary remarks: 


1, The region of stability in A, B, C parameter space must be symmetric with respect to the plane, 
A =B, since these parameters enter symmetrically into (24), 


= — 24, (B+ <0, 
and 
| 


2. If the signs of A and B are simultaneously changed, then the signs of all the roots of Equation (45) 
are reversed, Thus, if at point (A, B, C) to the real root, z; of Equation (45) there corresponded a real solution 
to System (24), then at point (— A,— B,— C), a real solution would not correspond to the root z;. In particular, 
bearing in mind the first remark above, we must conclude that if a real solution of System (24) corresponds to 
each real root, zj, of Equation (45) at the point (A, B, C) then there will be no real solutions at the point 
symmetric to this with respect to the plane A +B=0, 


_ 8, It follows from (30) that, for any real root, 24, of Equation (45), excluding the cases z=+ VC, we 
have 


4. Equation (45) has roots z= + VC, if and only if the corresponding point (A,B,C) lies on one of the 
parabolic cylinders 


B*=C. (47) 


If (47) holds, there are real roots, z =— A (z = — B), which follows from (30). 


5. If we change variables in (45), setting z= — — A, then Equation (45) with respect to will have the 
form: ¢ 


me + 252 (84*—C — 5) (449 4AC + B—A)+(AP—C) (4A) 


Since the constant term of this equation is positive for At C, then for all real roots of Equation (45) we 
have, for A** C, 


I] = @+4)>0 (48) 
t—1 


and for B*# C, by virtue of (46) 


k 
@ + B)>0. (49) 


where k is the number of real roots of Equation (45), By a simple transformation, Equation (45) may be put in 
the form 


For any real root, z, of Equation (45), there holds the inequality 


(A+B) +C+4+4AB>0. 


If A+B > 0, then it is not difficult to obtain from (51): 


At + 2(¢ +) 


From this last inequality, and taking (4€) into account, one may infer that for 


* The problem set here might be solved by taking as a necessary and sufficient condition the existence of at 
least one positive root of Equation (47a), However, using the corresponding criteria [3] leads to a really complex 
relationship which is unsuited for practical use. 


at) 


to 
rel 
(21 + A) (4, + B) >0. (46) pa 
als 
|| 
Aa 
qui; 
tha 
Aa 
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A+B>0, 2(C4 (52) 


to any real root of Equation (45) corresponds a real solution of System (24), But then we have, from our second 
remark, that for 


A+B<0, BY—2(C+7)>0 (53) 


(24) may not have real solutions, Conditions (52) and (53) determine, in parameter space, a region outside the 
paraboloid of revolution 


(54) 
also including its surface. 


It remains to consider the domain inside the paraboloid, characterized by the inequality 


<C ++ (55) 
In the sequel we shall always assume that Inequality (55) is satisfied, from which it follows that 


C42>0 
and, moreover, 


i 
C+ z+ AB>0. (57) 


Comparing the last two relationships with (38), we find that, at any point of Region (55), Equation (45) has 
at least one pair of real roots, 


In what follows we shall use a geometric presentation, We introduce the notation 


n=( (58) 
+ AB. (59) 


The investigation of Equation (45), together with (50) leads to the study of the intersections of the fourth 
order curve (58) with the straight line (59), Taking (56) into account, it is not difficult to find that the curve 


is tangent to the z axis for For 2 =0, yy attains maximum, ymax = OA = (2C + 1)* 


The graph of y, has the form shown in Fig, 1. Between the points A and Zp, and also between the points 
A and — zp, the curve has points of inflection for zing . 


When A, B, and C are in the Region (55), line yg cuts the y axis in the point M, above the origin, as re- 
quired by (57). Here, the distance from point M to point A on the y axis equals 


AM = +. AC B, (60) 


that is, it equals the constant term in Equation (45), Consequently, if C?— AB > 0, then the point M lies between 
A and B (line I- I on Fig. 1), Otherwise, point M, lies above point A (line I~ II on Fig. 1), 


It is also clear from (59) that the quantity A +B determines the angle between yg and the z axis, For 
A+B < 0, line yg is turned clockwise with respect to the z axis, and conversely for A+B> 0, 


Consideration of Fig, 1 leads to the following conclusions, 


| 
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1, For A+B > 0, Equation (45) has at least one positive root, Hence, in this case, System (24) has real 
solutions, since the positive root, z1, added to the positive quantity A (or B) gives z+ A > 0, but then by (46), 
20. 


2, If Conditions (38) and (39) are fulfilled and, in addition, AB <0, System (24) has real solutions, since 
there is also a positive root in this case, and one of the quantities,A and B, is nonnegative, It thus follows from 
(53) that for A+B <0, AB <0, Equation (45) can have only four real roots inside Region (55). 


c 


+ 


| 
yn | 

Fig. 1. 

Fig, 2, 


3. If A+B < 0, A® # C or B*3@ C and Equation (45) has two real roots, then neither of these roots leads 
to real solutions of System (24), Indeed, it is easily seen from Fig. 1 that at least one of the real roots, z;, of 
Equation (45) is negative for A+B < 0, Consequently, when added to the negative quantity A (or B), it gives 


% +A < 0, but then, by (48) and (49), for the second root also one will have zg+ A < 0 or zg+B< 0, Thus, it 
remains to consider the case when 


AB>0, A+B<0 (61) 


and the equation has four real roots, We shall consider the section of Region (55) by the plane C = const (Fig. 2). 
Here the curve yy is uniquely determined, Changing A and B within the curve formed by the section, we shall 
change the attitude of line yg and there rve the manner in which the roots of Equation (45) change, In 
view of (61) and the symmetry of the situation, it suffices to consider only the curve lying in the third octant 
below the line A =B, Suppose, for example, that C lies in the interval 1/4 < C <1, Here, the corresponding 
section (Fig, 3) will be intersected by the surface AB=C*, (It may be shown in general that for C > 
aw , Surface (54) does not intersect the 


surface AB =C?), 
fos s On Fig. 3 is also plotted the line 


A=—YVC, B=—YC, (62) 


obtained by the intersection of the plane C = const 


> 


Sd with the cylinder (47), Figure 4 gives the curve yy 
® corresponding to the given C, Here we limit our- 
selves to the right branch of the curve, inasmuch 
B80} x Why By) as we are only interested in positive roots, For given 
A = Ao, B =By (C is fixed), the line ys also takes a 


determinate attitude, It follows from Equation (59) 
that to a movement of the line 


A+ B=A,+B, (63) 


Fig. 3 in the plane C = const corresponds a displacement 
ig. 3. of the line yg parallel to itself along the y axis, but 
to the motion of the hyperbola 


AB = (64) 


corresponds a rotation of the line ys around the point M (z=0, y=C + Yq + AgBy). If the point (A, B) is moved 
along one of the lines (62) in the direction of increasing ||A |or |B | then the line yg turns in a clockwise sense 

around the point T (y = 4/4, z= ¥C), lying on the curve y;. This follows from the fact that, if one substitutes, 
for example, A = — VT in (59), one gets y»=(z— VC) (B — VC) + 4, i.e, an actual turning about the point 
T, That T lies on the curve y, is immediately verified. Thus, ifC = ¥,, the point T coincides with the point 


of inflection z;,¢ of curve yy, but if C > 1,, then T lies to the right of zinf, that is, in the region where 
<0. 


In order to establish whether or not a given positive root of Equation (45) provides real solutions of the 
system, it is obviously sufficient to establish whether z;, the given root, lies to the right or to the left of one 
of the lines 


t=—A,z=—B (65) 


We consider the attitude of the line ys (Fig. 4, line I— 1) corresponding to the point 1 (Ay— ¥C, C) 
(Fig. 3), At this point, AB < C® and, consequently, there are two negative and two positive roots, For the 
hegative roots, and zg, (zy + A) (z,+ A) > 0. 


For the root z,= VC = — B, we have z3+A > 0, since inside the sector under consideration, ~— A < — B, 


This means that for the second positive root, by virtue of (48) and (49), we have the inequalities 
M+A>0, 0 (zy 2%, since Ax B). 


Thus, at the point 1, both positive roots lead to real solutions, 


Continuously moving from point 1 to point 2 (Fig. 3) along line B = — VC, corresponding to a rotation of 
yg about the point T, we find that, at the intersection with the hyperbola AB =C?, the smaller, in absolute 
value, of the negative roots goes positive, but that the two other positive roots do not change their orientation 
with respect to lines (65), i.e,, as before they give real solutions to the system (Fig. 4, line I— II), Thus, at 
each point of the segment PK (Fig. 3), including the point K (A =B = — YC, C), we have real solutions for the 
system. Making use of the geometric representation, we may say that if we move from point P; (a, a) of the 
line A =B along the hyperbola AB= a* (a is some constant), intersecting the segment PK in the point Rj, we 
shall always have real solutions until we reach a certain point, Lj, coming after point Rj, but inside curve (54). 
The point L, lies on the discriminant surface, D =0, although at this point the System (24) still has a real 
solution, Along the remaining portion of the hyperbola leg under consideration there are no real solutions, For 
proof, consider, for example, point 2 (Fig. 3) and the line corresponding to it, I— Il (Fig. 4). If one moves 
from point 2 along the hyperbola (64) to line A =B, line I— II will rotate counter-clockwise about point My. 
During all this, there will be four real roots. In fact, if it is assumed that, in the process of moving, line I— II 
may occupy the position, shown in dotted lines on Fig, 4, then for the two real roots we would have in this 
position, we find that (z, + A) (zy + A) < 0, which contradicts (48), 


Consequently, the relative distribution of the roots of Equation (45) and lines (65) on the portion of the 
hyperbola PR; remains invariant, since it could change only with the intersection of line (62) or by the loss 
of two real roots, From this follows the existence of real roots in the portion PjRj, Moving from point Ry in 
the opposite direction along the hyperbola, one will always have four real roots and, consequently, real solutions, 
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until point Lj is reached, when line yg becomes tangent to curve y, (Fig. 3, striped line), At point 1j, obviously, 
D =0, and it lies within the curve since, as was proved above, for A+B < 0, System (24) does not have real 
solutions external to or on Surface (54), Thus, the required proof, We consider now a certain segment KN on 
the line A= — VC (Fig, 3). In this case, the root z3=+ YC = — A does not give real solutions, since B— A < 
<0. 


The negative root z; and the positive root z, which is lesser in absolute magnitude, give (zy + B) (zy + B)> 
> 0 (Fig. 5). For the remaining positive root z, we have z4+B < 0 and z,< %< 2, (Fig. 5). Therefore, on 
the segment KN, the System (24) does not have real solutions (excluding point K), 


” 


2=-8 


Fig. 5. 


Arguing as above, one may show that on any hyperbola (64) which intersects segment KN there are no 
real solutions, Moreover, it is possible to show that, if there is no solution at the point AgB,, then by moving 
along line (63) to the intersection with line A =B, still no solution will exist, Hence, it is not difficult to con- 
clude that, in general, for AB > C, real solutions of System (24) cannot be found in the sector under considera- 
tion, Consequently, the graph of the domain of stability has the form shown by the striped line in Fig. 3, and 
is determined by the relationships 


D=0, AB<C. (66) 


1 
The reasoning given above may be extended to cover all C > 4. For — 4 = C < 0, and for 0<C mar 
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the situation is somewhat different, since in the first case the plane C = const does not intersect the cylinder 
(47) and, in the second case, the point T (y = 4/4 , z = VC) lies on the portion of the curve yy for which yj > 0, 
i.e,, to the left of the point of inflection, 


It may be shown that all points of the strip, — + <Csx + , at which Equation (45) has four real roots, 


are found in Region (55), and at each of these points System (24) has real solutions, Thus, all possible sign 
combinations of the quantities A + B and AB have been treated, If A +B =0, then it is easily shown that a 
necessary and sufficient condition for the existence of real solutions to the system is the existence of real roots 
of Equation (45). 


Summarizing what has been said above, it is possible to assert that, for the existence of real solutions of 
System (24), it is necessary and sufficient that at least one of the following relationships not hold: 


A+B<0, AB>C>-+, (67) 
if Equation (45) has four real roots, and that the following relationship hold: 


A+B>0, (68) 
if Equation (45) has two real roots, 


In the final formulation, Condition (55) did not enter, by virtue of the remarks above, and also because 
for AB < C and Condition (65) satisfied, Equation (45) may have four real roots only in Region (55), Inequalities 
(15) and (16) are incorporated in the results obtained, entering into the stability criteria given above. 


It should be mentioned that by analogous means, stability criteria may be obtained for systems of arbitrary 
order [4], only if the problem may be tackled by quadratic equations of the form {I,4.6}. In truth, in the 
general case of an n'th-order system, the quantities A, B, and C, will not be rational functions of the coefficients 
of the initial System (1). 
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ON THE THEORY OF ONE-RELAY SYSTEMS 
V.S. Boiarinov and N, N, Leonov 


(Gor'kii ) 


In automatic control technology, frequent use is made of relay regulators of the “yes-no” type, 
Theoretical investigations of the dynamics of control systems supplied with such regulators have been made by 
a series of authors (see [1], and also (6, 7]). In particular, studies of some of the simplest models of such systems 
were carried out by Fifigge- Lotz and Klotter, using phase surfaces (space) [2-4], However, the results obtained 
in [2-4] are inadequate and, in a number of cases, erroneous,* 


In this paper we investigate one of the models of one-degree relay systems previously studied by Flfigge- 
Lotz and Klotter, here considering the model with various forms of relay characteristics, using line-to- line 


point transformations, To the study of this theoretical model, with the corresponding idealizations, is brought 
the theory of a series of second-order relay control systems, in which not only the coordinates, but also their 
first derivatives, play a part, Examples of such systems are the lateral stabilization of an aircraft autopilot 
where wing layout is V-shaped, a system of stabilizing autopilot longitudinal motion with a constant angle of 
attack for the jets, neglecting moment of inertia, etc, 


The system is described by a nonlinear differential equation 


with the relay characteristics shown in Fig. 1, Here, for generality, h is an arbitrary positive number (in [2-4], 
0<h< 1), The sign of the parameter k is arbitrary.* By means of Koenig-LeMaire diagrams a detailed 
investigation is made of the (two-dimensional) phase space of the given system and, in particular, new results, 
at least with respect to what has appeared in the literature, are obtained on the number of stable limiting 
cycles, There is also constructed a decomposition of parameter space into regions of qualitatively identical 
behavior of the dynamic system, Part 1 is due to N, N, Leonov, part 2 to V, S, Boiarinov, 


1, Case of Constant Spatial Delay 


Let the relay characteristics have the form A (Fig. 1), In this case the phase space of the system takes 
the form of two partially overlapping half-planes, joined along segments of the lines, y= § andy =~ 6 
(Fig. 2). If the representative point, describing the motion of the system, falls on these lines, the relay transfers, 
and the representative point moves over to the other half-plane, 


In half-plane I, the system motion is described by the differential equation 
y, y+ 2hy+a2=—1. (1) 


* For example, in works [2-4], for the dynamic system considered below, only one limiting cycle was allowed 
to exist, although there may be two of them, and sometimes even three, there lacked a decomposition of 
parameter space into regions of topologically similar phase spaces, etc. 

** A particular case of this problem was investigated, for § =0, by H, Bilharz (8), The case when k =0, but 
the characteristic of the nonlinear element has dead zones and loops, was treated by N, S, Gorskaia [6], 


Correspondingly, in half-plane I 


z=y, yt 2hy+2-=—1. (2) 
a) Leth <1, The solutions of the equations of motion, (1) and (2) will have the form: 
in half-plane I 
z=1-+e—" (acoswt + dsin wt); (3) 
in half-plane I 
(a, + sinwt), = (4) 


‘ The phase plot is symmetric with respect to the origin, since Equations (1) and (2) do not change when 
(x, y) goes over into (— x, — y), All phase trajectories may be divided into ordinary “plane” and for k < 0, 
“cylindrical” ones, The latter are found on the flattened surface during the generation of the cylinder which 
is generated by the two plane strips (between y = — § and y =6), joined in the region of the switching lines, 
¥ =—68,¥ = & . Movement along the cylindrical surface occurs if » changes sign at the moment the relay 
transfers, For ordinary ("plane") motion, % does not change sign at the moment of transfer, 


on 
Fig. 1. Fig, 2, 


The isolated points of the differential Equations (1) are (1,0) and (0,1), and are stable foci, These isolated 
points have finite zones of attraction, If the initial conditions lie in the zone of attraction of an isolated point, 
then the system will tend to a stable equilibrium state (x =1, y = 0, or x = — 1, y = 0), without transfer of the 
relay, An infinite phase space is unstable, since the equator of the Poincare sphere, in the given case, turns out 
to be a limiting cycle, 


Let the representative point begin moving (t = 0) from the point (xp, yo) on the line y=6 and, at time 
t=, land at the point (x,, ys) on the line y= — 6 (Fig. 2). 


Thena =2z—1, b = —[h(xo—1) + 


Setting t = r in (3) and in the derivative of (3), and substituting the values of a and b, we find the ex- 
pression for the corresponding function in parametric form (with r the parameter): 


-~8 
(1+ ( axe h) 


sin wt 


cak—2h+4—. 
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The time of transfer r varies within the interval (a A 1) for the “plane” case, and within the in- 


terval (0, r,) for the “cylindrical” case, ry and r, are determined by the condition that the phase trajectory 
and the line y= be tangent at the point (%, yy): 


from which 


The equations for ry, are of the form: 


(1 +8) (hsin — wcosat) = 0 (i = 4,2), 


a b 


Fig. 3. 


7 
S 


By virtue of the phase plot symmetry with respect to the origin, each time a transformation T* brings 
a point into the half-line L, (for which y = — 6), it is possible to consider instead the point symmetric to it, 
lying on the half-line Ly (for which y = 5), Thus the problem leads to the consideration of neon of 
the half-line Lg into itself by applications of the transformation T*, where the superscript bar indicatcs a transi- 
tion, after the application of transformation T*, to the symmetric point, Therefore, for the construction of the 
Koenig-LeMaire diagram, we investigate the curve x,= f(x»), described by Equation (5), For this, we find 


the expressions for and : 


dz, — o(1 +8) +(1 —8) (hsin wt + @ 008 wt) 

(8) 
sin%et [— 2h (1 — 8%) sin at +o (1 + 8)%e"* —a(i 

[o (1 — 8) + (1 + 8) (hsin at — @cos wt) (9) 


Investigation shows (see Appendix I) that the Koenig-LeMaire diagram will have five states (see Fig. 3). 
It is clear from Fig, 3 that, for k > 0, the following cases may occur: 


1) two limiting cycles (Fig. 3, d) (external stability and internal instability); 
2) one stable limiting cycle (Fig. 3, e) 
3) no limiting cycle (Fig. 3, c). 


For k < 0 (Fig. 3, a, b), there may be either one unstable limiting cycle, distributed on the flattened 
cylinder, or no limiting cycles, 


If k > 0 and the initial values lie in the zone of attraction of the isolated points, or internal to an un- 
stable limiting cycle (if such exist), then the system progresses to a stable equilibrium state at the positions 
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x=1 or x= -— 1, If the initial values lie in the region of attraction of a stable limiting cycle, then the system 
goes into self-induced oscillation, If k < 0, and there exists an unstable limiting cycle on the “cylinder,” then 
the representative point, leaving the cycle via some rotation about the “cylinder,” leaves the zone of attraction 
of the isolated points, 


The half- period of the periodic motion is determined from the equation * 


(h — k) sin ot — o3 cos ot = ash he + wd ch he. (10) 


For§—» Oand k <0, the “cylindrical” surface shrinks to the straight line y=0, The motion on the 
“cylindrical” surface, in the limit, goes over into a sliding regimen, 


b) Leth > 1, The investigation is completely 
analogous to that given above for the case h < 1. 
The isolated points of the system of differential 
Equations (1) and (2) are the vertices of stable junctions, 
The infinite phase space is unstable, since on the 
equator of the Poincare sphere there are four singular 
points, diametrically opposite in pairs, two of them 
unstable vertices, and two of them saddle points, 
Limiting cycles, besides unstable ones, do not exist 
on the surface of the flattened “cylinder.” In Fig. 4 
Fig. 4. Regions of topologically similar phase spaces is given a decomposition of k,h, § parameter space 
in the case of type A characteristics, I is the region into topologically similar regions of phase space 
with one unstable cylindrical limiting cycle, II is the (see Appendix II), Thus, in the case of a system 


region with one stable, and one unstable, “plane” with type A characteristics, depending on the para- 
limiting cycle, III is the region with one stable meters and on the initial conditions, the system 
“plane” limiting cycle, either will tend to the equilibrium states (1,0) or 


(— 1,0), or will develop autooscillation (Fig. 5). 


Fig. 5. 


Let the relay characteristic have the form B, As above, the phase surface consists of two half- planes, 
partially over- lapping, and joined along portions of the lines y=§ and y=-— 6. In half-plane I the system 
motion is described by the differential equations 


y+2hy+2=—1. (11) 
Correspondingly in half-plane I 

c=y, y+2hy+2=1. (12) 
Letting h < 1, we find expressions for the corresponding functions: 


1 —(i— 
(1 — 8) cos wt + (1+ (1 — 8) (h— hk) (13) 


* In this paper, ‘sh'= ‘sinh,’ ‘ch’ = ‘cosh'— Publisher, 
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from which 


— 8) + (1 + 8) (— hsin wt — cos wt) 

dz  —w(i+8)+ (1 —8)(-—-hsinat + acosatje™ (14) 
wt 2h (1 — sin wt + — 8)? + 

[— (1 +8) + (1 —8)(—/hsin wt + (15) 


As in the previous case, it is possible to have, not only the ordinary “plane” movement on the phase 
surface, but als » movement along a “cylindrical” surface, for certain parameter values and fork > 0, The 
isolated points of the system of differential equations (11)-and (12) are the stable foci type, with coordinates 
(1,0) and (— 1,0), These isolated points, for § < 1, do not lie in their region of definition. 


The time of transfer +r for movement of the representative point along the “cylindrical” surface, is 
within the interval (0, ry), where ry a, defined by the condition that the phase trajectory be tangent to 
the lines y=5, » =— 65. We distinguish three sub-cases, 


1, Tangency to y=6 at the point (Xo, yo) (Fig. 6, a), The condition gives % = 5 k (1+5)/c, and 
the equation for rj will be 


(wcos x, + hsinwt,) — — 4) = 0, (16) 
This condition should be used if § < 6* , where §* is determined from the equation 
—e 


2 ry for § =6*. Tangency occurs both with y =§ and with y= — 6, corresponding to the points 
(Xp, Yo) and Gry, yy). 
3, Tangency with » = —§ at the point (%, y,) which gives 


= 
The equation for r," has the form: 
—o(i+é)+(4 + 4) (— hsin wr, + w cos (17) 
and is used in case § > 6*. 
Yo) 
a 
Fig. 6, Fig. 7. 


The time of transfer r for “plane” motion changes continuously, with changing %», from 14" to , if 
6 < and from to ry, if > 6% where determined by the condition of tangency with y= — 6 at 
the point (%, ;). The equation for rgcoincides with Equation (17), but reer (Fig. 7), The time of trans- 


fer rz tends to = for motion of the point (x, y) along the line » = § out to infinity, 
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Investigations show (see Appendix III) that the LeMaire diagram may have the form given in Fig, 8, It 
is clear from Fig, 8 that, for k > 0, there may be the cases: one stable limiting cycle on the “cylinder” (Fig, 
8, b, c), one stable limiting cycle and two unstable on the “cylinder” (Fig. 8, d), one unstable limiting cycle 
on the “cylinder” (Fig. 8, f), one stable “plane” limiting cycle (Fig. 8, a), no limiting cycle at all (Fig, 8, g). 


Fig, 8, 


For k < 0, there may be the cases: one stable “plane” limiting cycle (Fig. 8,¢, j), one stable (external) 
and one unstable (internal) “plane” limiting cycle (Fig, 8, i), no limiting cycles (Fig, 8, h).* 


If the initial conditions lie in the zone of attraction of the stable limiting cycles, the system will develop 
auto-oscillation, If the initial conditions lie in the zone of attraction of the isolated points, internal to the 
unstable “plane” limiting cycles or parts of the "cylinder" not included between unstable “cylindrical” limiting 
cycles, then the system will move to stable equilibrium states at the positions x=1 orx=—1, For§—* 0, 
k > 0, the "cylindrical" surface shrinks to the line »=0, In the limit, motion on it tends to the sliding regimen, 
leading to the stable equilibrium state (0,0) to which, for § =0, the stable “cylindrical” limiting cycle degenerates, 
With decreasing §, the limiting cycles shrink, and the frequency of auto-oscillation increases without bound, so 
that it is possible to consider the equilibrium state (0,0), for § =0, as a degenerate limiting cycle, corresponding 
to auto- oscillation with an infinitely high frequency and an infinitesimally small amplitude, 


Assuming h > 1, we find the expressions for the corresponding functions: 


(1 + 8) ch — (1 — 8) e"* (1 + 8)(h— k) 


—(i —8)ch wt + (1 + (1 — 8) (h —k) 
¢ sh at — 3 + 


where w, = Vi? —1, c =k —2h+—+-, whence 


dey (1 — 8) — (1 + 8) ch + A sh (19) 
day (1 + 8) + (1 — 8) ch — sh wt) 

dz _ sh® — — 2h (1 — 8%) sh + (1 — — (1 + (20) 
de [— (1 + 8) + (4 — 8) ch — hgh 


* The cases a, d, f, g,h, andi of Fig, 8 are only possible for § > 1. 
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It suffices to consider the case §< 1 since for § > 1 limiting cycles do not exist, and the representative 
point, after no more than one transfer, leaves the stable positions x =+1 or x= — 1, 


2, 


Fig. 9. 


The time of transfer for “plane” motion varies within the interval (rq, oo) and, for motion on the “cylinder,” 
only possible for k > 0, within the limits from 0 to ry, The time ry is determined from the condition that the 
curve be tangent to the line » =§ at the point Xo, yo). From this condition: 


k(1 +8) 
— 


The equation for ry may be written as 
(1 + 4) (@, cha,t, + hshw,t,) — @, (14 — 8) 0, 


Investigation shows (see Appendix IV) that the LeMaire diagram can have the form shown in Fig. 9., 


From these diagrams it is clear that there exists one stable limiting cycle which, for some k , transfers 
from the “cylinder™ to the plane, or vice versa, 


Fer 6 <1, any initial values lead to a stable limiting cycle and, for § > 1, to stable equilibrium positions 
(vertices) with coordinates (— 1, 0) or @1, 0), 


For §— 0, the “cylinder” generates a sliding regimen, and the stable “cylindrical” limiting cycles tend 
to the stable isolated point (0,0), 


Figure 10 gives the decomposition of parameter space into regions of topologically similar phase spaces 
(see Appendix V), 


a / 
? 
| 
net 
55 
= = ¢ 
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Fig. 10, Regions of topologically similar phase spaces for type B 
characteristics (Fig, 1), 1) The region with one stable “plane” 

limiting cycle; 2) the region with one stable and one unstable " 

“plane® limiting cycle; 3) the region with one unstable “cy- 

lindrical” cycle; 4) the region with one stable and two unstable 
“cylindrical® limiting cycles; and 5) the region with one stable 
"cylindrical® limiting cycle (equations for the surfaces S;,..., Be 
Sg are given in Appendix V), _ 
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From all that has been said above, it follows that if a relay has type B characteristics (Fig, 1), the re- 
presentative point, for 6 < 1, inevitably tends to a stable limiting cycle, In all other cases (type A char- 
acteristics, and type B for 6 > 1), despite the presence of stable limiting cycles, the representative point, with 
corresponding excitations, is not guaranteed to fall within the zone of attraction of an isolated point, 


In conclusion, we summarize the expressions from which the parameters for auto- oscillation are deter- 
mined, 


Type A Relay Characteristics (Fig, 1), The period of auto-oscillation, T = 27*, where r* is determined 
from Equation (10), Knowing r* , we can find x, from the first of Equations (5) after which, determining a and 
b from Equation (3), we find the system motion with auto- oscillation, 


Type B Relay Characteristics, The period of auto-oscillation is T = 27r*, where the half-period r* is 
determined from the equation 


(hk — k) sin ox + 8 cos wt = w sh ht — w ch he. 
Knowing r*, we find x from the first of Equations (13) after which, determining a and b (a =x, +1, 
b = (% +1) + / @),we find the system motion with auto- oscillation: 
z=—1+e-"(acoswt + bsin wl). 


2. Relay Characteristics With Zones of Instability 


Let the characteristic of the nonlinear element have the form of C (Fig, 1), The problem is described by 
a system of differential equations of the form 


(9), 
1 for >4, 
0 for (21 
—1 for po 


As above, the solution is found by considering certain linear differential equations with successive “aligned” 
stationary integrations, starting from the requirement that x and y be continuous at the points of discontinuity of 
the right side of (21), In the problem considered, such an alignment leads to the lines y = 4 6. 


The phase space appears as a plane, decomposed into three 
regions, Gy, G, and Gg, joined at the lines y=§ and y= — 6 (Fig. 11). 
Since Equation (21) is invariant when (x, y) goes over into (— x, — y), 


\ 
S fs Ce F the phase space is symmetric with respect to the origin, For § < 1, 
Os there are three equilibrium states at the points (— 1,0), (0,0), and (1,0), 
of the stable vertex or focus type, encircled by their respective zones of 


attraction, For § > 1, there is only one equilibrium state, (0,0), with 
Fig, 11 one zone of attraction, The respective point, starting its motion in 

i? region Gz will either asymptotically approximate to the equilibrium 

position (0,0) if the initial state of the representative point lies in the 

zone of attraction of the point (0,0) or, after a finite time, it will fall on one of the lines y= +6 and trans- 
ferred by it, will move along the trajectory adjoining this line on the other side, following the corresponding 
other differential equation, The phase trajectory of region G, is tangent to the line y = § in the point with the 
ordinate 


The phase trajectory of G, is tangent to this line at the point with the ordinate 
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45) 


Therefore, the aligning initial conditions are impossible along segments of the line » =6§ for values of y in 
the intervals a 
3)k 
> 


since, on these segments of the line y = 6, the trajectories adjoining the line are directed in different directions 
(similarly for the line y = — 6), In these cases there occurs a sliding motion of the representative point along 


the lines satisfying the equation %= +x + 6). The sliding regimen is stable for k < 0 and unstable for 


k > 0, If the sliding regimen is stable, the representative point moves along the lines y = 4 § to the boundary 
of the zone of attraction, then enters the zone of attraction, If the sliding regimen is unstable then, withan 
arbitrarily small excitation, the representative point leaves the lines » = 4 6. 


Leth < 1, Then the solution of differential equation System (21) has the form: 


(a, cos wt + sin wt) for 


xz = e— (a, cos wt + b, sin wt) for <3, (22) 
+ bssinwt) for — 8, 


Initially, at time t =0, let the representative point be at point (xp, yo) of the half-line Ly and after time 
Ty let it fall on point (x;, y,) of half-line Ly via trajectories in region G, (transformation S*), Then 


a, =8—1—kyp, Wt B— —h 


The corresponding function for the transformation S* (in parametric form) is written as follows: 


1—8/ coswt,; —e"™ 
Ch ( o+h—k), 


sin wt, 


where cy =1 — 2hk+k’, 


The time of transfer ry varies within the interval if 6 < 1, and in the interval (0, ) 
if § > 1, The equation for ry is determined by the condition that the phase trajectory touch the line » = 6: 


= — hsin — wt, + we = 0, (24) 


dy, —@+(hsin ot, +o cos 


fo (25) 
dyo + (hsin wt, — @ cos @t;) hh 
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3k 
whence 
In investigating the corresponding function, we find 
| 


sin? wt; 2 (hsin ot, — sh ht) — %siniwt, fs 


In the interval (0, Fj), the following inequalities hold: 


fe<Oand <0, >0. 


The curve y;(yg) has the asymptote 


(2h—k)(ite ®), 


Now let the representative point leave the point (x, y,) of the half-line Ly at time t =0 and, at time 
t= Tq, fall on the point (xs, ys) of the half-line L, via trajectories in region Gg (transformation E*), With this, 


§—ky, b, = nt . The corresponding function has the form: 


na + (= 


sin 


— COS — 


sin 


n=+2( w+h—k). 


(27) 


The time of transfer r, varies in the interval (0, where Tz Since the equation 
1 
for T; has the form: 


— hsin + wcos wt, + = 0. (28) 
It follows from Equation (27) that 
dy. @ + (h sin wT, + @ cos wT.) 
w+ fs ’ (29) 
sin? wt. 2(hsin wt, + w sh sin?wt, fe (30) 


In the interval (0, 73), the following inequalities are satisfied: {,> 0, f, > 0, and fg > 0, Consequently, 
and 
dy 


Further, let a transformation S take points of the half-line L into points of the half-line Ls ( » = — 5) 
along trajectories in region Gs, and a transformation E~ transfer points of the half-line Ls into points of the 
half-line Ly (y = §) along trajectories in the region Gp, We assume that the points (xo, yo) (x4, ys), (Xa yo), 

» lie outside the segments of sliding motion, The traisformation T=S*E'S“ takes the half-line into 
itself, 


By virtue of the symmetry of the phase plot with respect to the origin, it is always possible, when a point 
is brought by the transformation T =S*E* to a point on the half-line L», to consider instead the point symmetric 
to it, lying on Lo, Therefore, the problem leads to the consideration of transformations of the half-line L» into 
itself via applications of transformations, T = S*E*, where the superscript bar indicates that there is a transfer to 
the symmetric point after application of the transformation S*E*, 


If the transformation T is not applicable to certain points of the half-line Lo, this means that the 
corresponding representative points lie either in a zone of attraction or on lines of the sliding regimen, 


Figure 12 gives some of the LeMaire diagrams constructed in accord with Equations (23) and (27), It 
is clear from Fig. 12 that, for k < 0, limiting cycles are not possible, For k > 0 there may occur either stability 


> 0, The curve y2(y,) has the asymptote nent (2h— k). 
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in the large, or two limiting cycles (stable inside, unstable outside) or one stable limiting cycle, If the initial 
motion lies within the region of attraction of a stable limiting cycle, the system will fall into auto-oscillation, 


4 
% 
(8<1;2h-x<0) 


Fig. 12, 


Its period will be T = 2 (ry + rg) where ry and Tr, are the times of transfer corresponding to the fixed points of 
the transformation, If the initial values lie in the zone of attraction of an isolated point or exterior to an un- 
stable limiting cycle (if such exist), the system goes to a stable equilibrium state at the position (1,0), or at 
(— 1,0) or, finally, at (0,0), 


We now find a decomposition of k,h, § parameter space into regions of qualitatively identical dynamic 
system behavior, In accord with Fig, 12, the condition that there appear a semi-stable cycle (the condition of 
mutually tangent curves) determines the boundary dividing the region in which two limiting cycles exist from 
the region with no limiting cycles, These conditions may be written as 


or, in expanded form, 


cos wt, — +h—k cos wt, 


(1 — h~k de, 


sin wt, sin wt, (32) 
— + (hsin wt; + cos e~"™ + (—hsin wt, + cos wty) 
+ (h sin at, — w cos wt) e"** + (h sin wt, + @ cos wt) (33) 


Equations (31) = (33) are the parametric equations for the limits sought (with ry and r3 as parameters), 


From the condition for the disappearance of unstable cycles (one of the points of intersection of the curves 
lies on the boundary of the zone of attraction) we determine the boundary dividing the region in which two 
limiting cycles exist from the region in which one limiting cycle exists, These conditions are written as: 


for§<1 
Yo (™2) = — Yo Yi = (1); 
for §> 1 
or, in expanded form: 
wt, — cos wt, + we" = 0, 
cos wt, — Cos wt, + 
(1—4)@ = OTs (34) 
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k Bh cos wT, + 


sin 


for > 1 


— hsin wz, + w cos wt, + = 0, 


(i — cos wt, + — 008 wt, — 


sin SiN 


(35) 


k= h(1—2) + 


Sif WT; 
In the expressions given above it is necessary, forh > 1, to replace w by iwy, where w= vh — 1, The 
decomposition of parameter space for the case under consideration is given in Fig, 13 (see Appendix VI), 


If the relay has type D characteristics (Fig, 1), then the problem is described by the following system of 
differential equations: 


y, 
y+ 2hy +2=F 
—1 for p>4, 
(36) 
1 for p< — 3. 


The investigation of this problem also leads to point transformation of a line into itself, By virtue of the 
symmetry of phase space it is possible, here as previously, to limit oneself to half the point transformations of 
a half-line Lo into the half-line L, (transformations T, =S} Ej’), If we let b=1 —6, and by = — (1+ 6), it is 
then easy to show that the transformation T, = S{Ej‘ is obtained from the transformation T =S*E*, in the case 
of type C characteristics, by replacing b by by (for b < 0) in the corresponding expressions, 


h 
2 
— 7 
Fig. 13, Regions of topologically similar phase Fig. 14, Regions of topologically similar 
spaces in the case of type C characteristics (Fig. 1). phase spaces in the case of type D character- 
1) The domain with one stable and one unstable istics (Fig. 1), 1) The region with one stable 
limiting cycle; 2) the region with one stable limiting cycle; 2) the region with one stable 
limiting cycle, and one unstable limiting cycle, 


Independently, of the magnitude of §, there exists only one equilibrium state, (0,0), with its own zone of 
attraction, since the isolated points, (— 1,0) and (1,0), do not lie in their domain of definition, For k > 0, the 
system is stable in the large, For k < 0, limiting cycles may exist, The decomposition of parameter space is 
given in Fig, 14, 


For type C characteristics, the parameters of auto-oscillation are determined from the following ex- 
pressions, The period of auto-oscillation T=2(rf + 7%), where ry* and 1,* are defined by the system of 
Equations (31) and (32), Knowing ry and 12 we find yp and yy from the first equations of (23) and (27) after 
which, determining a, and b, (also ag and bg) from (22'), we find the motion of the system under auto-oscillation 
from Equation (22), In the case of relays with type D characteristics, the parameters of auto-oscillation are 
also found without difficulty, . 
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APPENDIX I 


Investigation of the Function Corresponding to (5) 


In the investigation of Function (5) we first of all determine the sign of = and fam - Since sinwrt 


is positive in the interval (0, r,) and negative in the interval a, 74), it is easily verified that the numerator 
of Expression (8), — «(1 + —8)(hsinwt + wcoswt)e , is negative in the interval (0, =) and 
that the denominator of Expression (8), «(1 — 8) + (1 + 8)(hsin at — wcos wt) eis negative in the interval 
(0, r2) and positive in the interval( =~ ° ry. Further, it is easily seen that the expression within the square 
brackets in the numerator of Formula (9), — 2h (1 — 8%) sin wt + @(1 + 8)*e"* —o(i—)*e—"", is positive 


a ) since wsinh r > hsinh wr. 


in the interval (, 


It thus follows that > in r) and <0 and thatthe sign of in the 


intervai (0, ) depends on the sign of c, i.e,, in the given case, on the sign of k: a > 0 for k< 0 


and > <0 fork >0, We then find that, for the “plane™ case, the curve x = x; (xq) has the asymptote 


(or 


mh rh mh 


For the “cylindrical” case, the curve xX, = x(x») has the asymptote (for r—> 0) xy =X9 — 26 ke. 
We also note that 


= (Te) == 6D. 


All the above assertions determine the form and orientation of the curve x, (x9), 


We mention in conclusion that the equation for the half- period, (10), has no solution in the interval (=, 
m1) . Consequently, on the corresponding Koenig-LeMaire diagram the line x;(x9) will not ‘ntersect the line 


X= — X (for < T< Ty), 


APPENDIX II 


Decomposition of k, h, 6 Parameter Space 


The boundaries of the topologically similar regions of phase space, forh < 1, are determined by the 
following conditions (see Fig. 3), 


1, For k < 0, from the condition that unstable limiting cycles appear, This condition determines the 
boundary of the surface dividing the region where there are no limiting cycles from the region where there 
exists one unstable limiting “cylindrical” cycle, The corresponding equations are written as: X,= — X= % 
(X, is determined from Formula (6)]) and, in expanded form, 
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w (1 —8)e—"* = (w cos ot — hsin wt) (4 + 8), 
(h — k) sin ot — 8 cos wT = w sh At + @8 ch ht. 
2, For k > 0, from the condition that unstable limiting cycles appear (from the condition that the tangent 
to the curve x, (x4) on Fig, 3 be the bisector of the right angle), This condition determines the surface dividing 
the region in which no limiting cycles exist from the region in which two cycles exist, The corresponding 


equations take the form: 


— ty 


or, in expanded form, 
+ (1 + 8) (w cos wt — hsin wt) = (1 — 8) (h sin wt + cos wt), 


(h — k) sin at — w8 cos wt = w (sh At + 3 ch ht). 


3, Fork > 0, from the condition that the unstable limiting cycle join the boundary of the zone of attrac- 
tion of the isolated points (— 1,0), (1,0), This condition defines the surface dividing the region in which there 
exist two limiting cycles from the region in which there exists one limiting cycle, The corresponding équations 
coincide with the equations for k < 0: x4 = Xo =X; 


For h > 1, there exists only the continuation of the surface I, Its equations are 


©; (1 — 8) = ch — h sh (1 + 8), 
(h — k) sh wt — «8 ch @yt = @ (sh ht + 3 ch hr), 


where = Vie—T. 


APPENDIX III 


Investigation of the Function Corresponding to (13) 


For investigating Function (13), we determine the signs of “{- ana SY, As was stated earlier, the 


time of transfer r varies within the following intervals (see table below) where 1," and ry" are determined 
from Formulas (16) and (17) : 


3 2 Motion on "Plane" It is easily verified that the numerator of Expression 
& < §*, is positive in the interval (x: =) and negative. 
@ in the interval (0, ry) but, for § >6*, is negative in both 
0. =, Te the intervals (0, 7") ana 


As for the denominator of Expression (14), — @ (1 + 
+8) + (1 — 8) (@cos —hsin wt) e"* it is negative for 
& < 6° in the intervals (0, r) and (rs - +} » negative for §> § in the interval (0, ry") and positive for 


& >6* in the interval r). 


We now introduce the notation f(r) for the expression in the square brackets in the numerator of Formula 
(15): f(t) = — 2h (1 — 84) sin wt + (1 — — w(1 + 


We wish to determine the sign of the derivative, /’ (t) = hw { — 2 (1 — 8) cos wt + (1 — 8%) e"* + (1+ gyre}. 
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We immediately convince ourselves that, for § > 1, coswr> 0 and for § < 1,coswr< 0, the deriva- 
tive > 0, 


We proceed as follows in order to determine the sign of f(r) when § > 1, coswr < 0 and when § < 1, 
cos wt> 0, Let us assume that, in these cases, f* (r)> 0, Then 


2 (1 — 8%) cos wt < (1 — 3%) 4 (14 


Here, the left member of the inequality is positive, It is therefore permissible, without changing the sense 
of the inequality, to square both members, leading to the expression 


[(1 — (4 + 4 4 (1 > 0, 


the correctness of which is obvious, If now f'(r) <0, the sign of the preceding expression would have to be 
reversed, which is impossible, Thus, f'(r) > 0 for all r, In the intervals (0, 7y") and (rq',#/ w) there exists 
ar =r", such that f(r") =0, It is easily verified that f(r) > 0 in the imervais 74"), and (r’, 


+) and f(r) <0 in the intervals (0, rj), (0, r*) and (ry', r"). It follows from what has been said that 


0 in the intervals (0, 14") and (0, 74"), and < 0 in the intervals [r,', and 


With the condition that § < §* in the intervals (0, r4') and (r,', r'), or for § >§* in the interval (0, r‘) 
d*zx, 1 
0, if k>O 0, k<0. 
> anf < if k< 


With the condition that § <* in the interval or for § > §* inthe intervals (r,', 7") 


and Can > ok 0, and < 0ifk> 0, 


Incidentally, we may mention that 


The curve x, = xX, (xq) has the asymptote (for ry—> 0) 
28 


and the asymptote (for r—> +) 


2h —k 


It is necessary to consider in greater detail the 
case k> 0, & > 6* (Fig. 8). In this case, a variation 
a b i in the system parameters may occasion a reversal of 
a the stability of the limiting cycles, In order to follow 
a this out in detail, it is necessary to consider not only 
the point transformation Xq = x4 (%), but also the in- 
Fig, 15, verse transformation, xq = (x4), 


In this case, the corresponding function has the 
form: 


— cos wt — (1 +8) 
csin@t 
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dz, dz dx 

dz» (t,) = 0, = 00, =o. 
| 
|| 


— — (1 — 8) e™ 


esin @T 


As a result, the diagram shows the curve xj = xj (x) as symmetric to the curve x = xq (x9) with respect 
to the line xy = — X%, The different cases which might arise here are shown in Fig. 15 (cases a, b, c), 


It is necessary to show that the curves of Fig, 15 are not intersected in more than three points, Additional 
points of intersection, over and above those shown in Fig. 15, may occur in sections where both curves simul- 
taneously are convex upward, Consequently, it is sufficient to consider the interval (r», rT) of variation of r, 
in which the second derivatives of x, (x9) and xj (x) with respect to X» are negative, ro is the value of + for 
which 


and T, as above,is determined from the condition of contact: 
— + 8) + (1 — 8)(— A sinwt + wt) 


Wy ta) — If x (xp) = 0 at not more than three points, and con- 
sequently, if St =0 at not more than two points, then the curves in Fig. 11 ate intersected in not more than 


three points, We now show that the curves fy) = and f(x) = intersect no more than twice in 


(ro, 7). Such will be the case if, in (ro, T), > 0 and 0. 


We show that 


a? 
> 0 in (to, 2). 
We have 
dt 
dz dt \ dz? ) dz, 
But, since 


_—o(1 +8) + (1 —8)(—h sin wt + wr) 
dt a sin*wt 


<0 


in (0, F ), it is necessary that the expression 2(#y) be less than zero in (ro, T ). 


We introduce the notation 


— + 8) +- (1 — 8) (— Asin wt cos wt) e"* 
sin wt 


= fs(t), 
(1 — 8) h sin wt + (1 — 8)? — (1 4-8)? = fy (4): 


With this, 


a ry . 
—fe (3) = 13 (7), 
dfy(t) d sd?x d*z, 
Here, +> 0 forall r, since {,(r) coincides with f(r) (see this Appendix, above) fy < 0 in (0, T) 


118 


= 


a 


and < 0 in (ro, TF). 


We consider the derivative 
diy + 8) cos wt — (1 — 8) (sin*et— hasinet coswt+*cos*wt) 


dt sin’wrt 
With this, 0 in(o, and in par ofthe interval 4). dr 
may be less than zero, From > et follows 
14-8 sintwt — hw sin@t cos wt + 
1—8 w? COS @T 
Substituting the expression for —- in a ,» we obtain 
Is 
dz? >0, i. e., in (To, 
d(d 
because in (ro, T ), a” <0, Consequently, from (A) 0 in (ro, T). 
We now show that o <0 in(ro, 7). 
We have 
d*f, d d*z,’ dt 
dz? = ae( dz? 
But, since 
dz dz, —w(i—8)+(1+8)(hsinat+ wt)e™ 
a” asin? wt 


in(o, 116 ta he exprenion be less than zero in (ro, rT). We introduce the 


notation 


= 8) +4 +8) (sin wt + w cos wr) 
sin wt 


Then 
a 
— — ff, (7) f5(*), 


Ga) 
<0in (%, 7). 


But, a ( *\, and 
0 for allt, fe>0in(0, =) 


We consider the expression 
w? (1 — 8) cos wt — (1 8) e~** + hasinwtcoswt + cos 


dfs _ 
dt sin*’wr 


d 
Here <0 in (0, 7*), where “>>. 
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From = 0 it follows that 


1—8 sin*wt-++hosin wt coset + w*cos*wt 
COS WT 


Substituting the expression for ak in Sat we get . 0, whence r*> T, i.e., in(ro, T), 


0 and, from (8), 


In case g (Fig. 15) there is a stable limiting cycle, At the moment of transfer, at point O, via the curve 
xq = xq (x) and at point O, via x, = x, (xq), there appear two unstable cycles, between which are found stable 
ones, Further, with variation of the system parameters, these unstable cycles approximate to the stable ones 
and fuse with them, As a result, there appears but one unstable limiting cycle, 


APPENDIX IV 


Investigation of the Function Corresponding to (18) 


For the investigation of the function corresponding to (18), we determine the signs of —% and ae" 
We designate by f,(1r) the numerator of Expression (19), /; (t) = «1 (1 — 8) — (1 + 8) cha;t + h sh 


and thus find that fj (r)> 0 for any positive r, Let f,(ry)=0 for r= 7, Then fy < 0 for 0 < r < 14, and 
fy > 0 for ry < r < @, 


Let us designate by f, (7) the denominator of Expression (19), It is easily seen that f,(0) < 0, fg (r) < 0 
and, consequently, f, (r) < 0 for0 < r < 


Further, in Formula (20) let 


f(t) =— 2h (1 — 8%) sh + (1 — 8)% 
We then find that 


f(m)<0, 


and, moreover, =0 for r satisfying the relationships and ae +o)*" where 
Since f(t”) =(1 — 8)%h + (1— )< 0, f(r) consequently changes sign only once, for some 
r’> T,, i.e., in the case of “plane” motion, From this, it is not difficult to see that 

for 0<t< 7,. <9 for <tT< =, 


d?z, for O<k<h—w, or k>h+,ininterval (0, 7), 7), 


or or k>h+, in interval 
dr? for k<0 or jninterval (0,71) 7°). 


The curve x4(x) has (for r—> the asymptote 


h—k— — 4814" 


c 
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y_ 
a (Ga )<° 
0 
dz k<0 Or h—w<h<h+a, in interval (70, co), 
| 


For r—> 0 the asymptote is 
28 


APPENDIX V 


Decomposition of Parameter Space 


For h < 1, the conditions defining the boundaries of the topologically similar regions of phase space will 
be the following (Fig. 8). 


1, The condition of transfer of a stable limiting cycle from the “cylinder” to the “plane” or conversely, 
gives thé equation for the boundary of the surface S;: xy = — xX) = — % (here, Xp= § — kus) or, in ex- 
panded (parametric) form 


(1 8) (w cos wt + h sin wt) = (1 — 8) 
(1 — 8) cos wt + (8h —h + 2k)sinat (wr <n). 


2, Thc condition that two unstable limiting cycles appear on the “cylinder” gives the equation for the 


surface Sg: x = xy =%j, Xo (here, % =— 6 or, in expanded form: 


(1 +8) = (1 — 8)(— hsinwr, + cos 
(a1 
(1 +8) cos wt, + (h + 8h — 2k) sin wt, = (1 — 8) wel, 


ay hv —ht 


3, The condition that three cycles on the “cylinder* merge into one unstable one gives the equation 


for surface Sy: x4 = — Xp, =1 or, in expanded form, 


280 cos wt + 2(h — k) sin wt= w — B)e"* —(1 + Sx) 


2a — (1 + 8) (h sinwt + w cos'wt) e~"* = (1 — 8) (—h sin wt + wcos wt) e”*. 


4, The condition for the disappearance of an unstable limiting cycle from the “cylinder® gives the 
equation for surface Sg: x, = — = % or, in expanded form, 


(4 —8)(—Asin wt + wcos@t) = (1+8) 
(4 + 8) cos wt + (8h + h — 2k) sin wr = (i — 8) 


5, The condition for the emergence of an unstable “plane” limiting cycle gives the equation for surface 


Sg: %, or, in expanded form, 


280 cos wt +2 (h—k) sinwat =o [(i—8) wt < 


208 + (1+ 8)(hsinat + wcoswt)e—"* = (1 — 8) (@ cos ot —h sin wt) 


6. The condition that an unstable "plane" cycle degenerate to the boundary of the zone of attraction of 
the isolated points (— 1,0) and (1,0) gives the equation for surface Sg: x, = — Xj = %, coinciding with the 
equation for surface S 4, 
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In the equations for surface Sg, wr varies in the interval (wm < wr < 2), Forh> 1, there exist (Fig, 9) 
only the extensions of surface S, 


(1 + 8) (ach sh ct) = (1 — 8) 
@; (1 — 8) ch et + (8h — h + 2k) sh wt = (1 + 8) aye" 


and the plane § = 1 divides the region devoid of limiting cycles from the region in which there exists one 
limiting cycle ("plane" or “cylindrical"), 


APPENDIX VI 


Construction of Parameter Space Decomposition 


To construct the decomposition of parameter space into regions of qualitatively identical dynamic system 
behavior, we subtract Equation (32) from Equation (31), The result is 


1—8 cost, —chht; __ 008 wT + ch hte (a) 
sin SiN wT, 
Up |p u, |, Equations (33) and (a) may be written in the form 


ua (h, 13) = (hia), (b) 


t, (h, T2) = (h, 73). (c) 


The limits of variation of ry and r, were stated above, We 
fix h, and construct the curves corresponding to the left and right 


members of Equations (b) and (c) for various a = apd (Fig. 16). 


Let the curves ug(r3) and v2(r3) be plotted on transparent paper, Superimposing the plane ug(vq), rz on the 
plane uy (vy), Ty, we slide the planes along the axis of abscissas seeking the points of intersection of the con- 
structed curves, The system of Equations (b) and(c) will have solutions rj and rif the points of inter- 
section of the curves ug with aw, and v2 with v, lie on one vertical, Determining rj] and r} from Equation 
(a), we find the points which lie on the surface dividing, in 5, k, h space, the region in which two limiting 
cycles exist from the region where no limiting cycles exist, 


For finding the boundaries dividing the region in which two cycles exist from the region in which one 
cycle exists, we take Equation (34), with h fixed, in the form 


Ue (T2) = au, 
and, correspondingly, Equation (35) in the form 


Ue = (7). 


Then, in a manner analogous to that given above, we find the points of the bounding surfaces. 


Fig, 16, 
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DETERMINATION OF PARAMETERS FOR CORRECTIVE DEVICES IN 
LINEAR SERVOSYSTEMS USING GIVEN GENERALIZED PARAMETERS 


M, M, Kreimerman 


(Moscow) 


An analytical method is described for determination of parameters for correcting 
devices in a linear servosystem when the values of some of the generalized parameters 
of the system are given, such as: the asymptotic crossover frequency, the phase margin 
of stability and the relative quality factor, A table of basic formulas for calculations 
are given, and examples of computations are carried out for series and parallel correct- 
ing devices and also for their combinations, 


INTRODUCTION 


In the present paper a procedure is given for selection and determination of parameters for correcting 
devices, which makes it possible to solve the second stage of the synthesis problem for servosystems as formulated, 
for instance, in[1). The problem treated in the present paper, however, differs from the problem of [1); here it 
is.not an over-all optimal characteristic of the entire system which is given, but only optimal values of some of 
the most important generalized parameters which affect the precision of the system, In servosystems, which have 
typical logarithmic amplitude characteristics (1, a, c.), such parameters a the asymptotic crossover frequency 


We, the phase margin of stability ¢,, and the relative quality factor D, = > (Dis the quality factor, or 


velocity amplification factor of an astatic system of the first order), Thus, Pe. problem consists in determining 
the transfer function and the parameters of a correcting arrangement, when the following are given: the shape 
of the typical 1.a.c., the transfer function of a certain part in the block diagram of the servosystem whose 
components have been selected for technical reasons, and the numerical values of w., @, and Dg. Let us 
emphasize that we, %, and D, in turn are selected in such a way that at the given magnitudes and character 
of the input pulse, the over-all error of the system should be at a minimum, 


When selecting one of the four types of l.a,c, [1] it is necessary to take into account the following: 
1, Types I and Il of the l,a.c, are the simplest for practical realization; 


2, The error in a system having a 1l.a.c. of type II or type IV can be made smaller in magnitude than 
for a system with a l,a.c, of type I or type II; 

3, In the presence of perturbation input only, the smallest error can be obtained with a system having a 
l.a.c, of type 1; 


4, If the system is under the influence of both controlling signals and perturbations, and the relative 
quality factor is large, with l.a.c, of types Il and IV a somewhat smaller minimal error can be obtained 
than with types I and III, 


The block diagram of a servosystem with series and parallel correcting devices usually has the form shown 
in Fig, 1. 


1% 


The transfer function of the open circuit loop for a servosystem of the above-described kind, as is known, is 


K, (p) (p) (1) 


where K,(p) = Ky(p)Kg(p) is the transfer function of the open circuit loop of the given (uncorrected) part of 
the system, K,(p) is the transfer function of the series corrective device, Koc(p) =Kg(p)K'k(p) is the open 
circuit transfer function of the local feedback loop, K*,(p) is the transfer function of the parallel correcting 
device, K,(p) is the transfer function of the system links not encompassed by the local feedback loop and 
K,(p) the transfer function of the system links included in the local feedback loop, 


e(t) 


A 


Kip) 


Fig. 1. 


It is necessary to select K(p) and K (p) in such a way that the over-all open circuit transfer function 
of the system 


K (p) = K; (p), (2) 
where we, @c, and D, have the chosen optimal values, 


Here K,(p) is an open circuit transfer function for the system, which has a 1.a.c, of either typeI, II, II, 
or IV of those given in paper [2]; thus i =I, II, Il, or IV, 


In special cases 
K (p) = Ky (p) K, 
when only series correcting devices are used in the system, and 
4) 


| = TF 
when only parallel correcting devices are used, 


1, Basic Relationships for Computations 


To find the basic relationships needed for this computation let us investigate the region of each of the 
typical l.a,c, [2] in the neighborhood of the asymptotic crossover frequency we. This will permit us to find 
the connection between the parameters which are included directly in the open circuit transfer function and 


the generalized parameters (w,, ¢c, and D,). Approximation equations connecting the two sets of parameters 
are as follows: 


for a l,a,c, of type I: 


for a l.a.c, of type I: idl 
22 (1—y/ Se) 4 (5) 
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for a l,a,c, of type Ill: 


for a 1,a,c, of type IV: 


Here 6 (g) is a linear function of g, with a negative slope: 


Be) 


In Table 1 there are given approximated numerical values of this function for some values of the phase 
margin of stability 


TABLE 1 


%c» deg. 20 2% 30 35 40 45 50 55 | 


B(@) 0.60 | 0.56 | 0.52 | 0.48 | 0.44 | 0.39} 0.35} 0.3 | 0.28 


Deduction of Equations (5) is given in an appendix to this paper. 


For a servosystem with large relative quality factor | V 2 > 1) Equations (5) are simplified to: 


B+ (9). 


To solve equations of the type of (5) or (7), it is necessary to have an additional condition given, This 
conditic .. can be, for instance, that the ratio ws3/we2 for a given ¢_ shall have a certain minimum value, or 
in other words, that the frequency band width of the l.a.c, with a slope of 6 db per octave (or 20 db per decade) 
should have a given minimum value in the neighborhood of the asymptotic crossover frequency @¢. In this 
case the servosystem will have a slowly changing phase- frequency characteristic of the open circuit loop near 
the crossover frequency (Fig. 2), which will assure the constancy of y., one of the basic parameters of the system 
when the amplification factor of the system is changed, The latter is important from the point of view of 
maintenance needed during operation of the system, In addition, the fulfillment of this condition is also of 
practical importance in tle following respect: when correcting the characteristics of the servosystem, for in- 
stance, by means of derivative loops, the quantity ws/wz is numerically equal to the loss in the quality factor 
of the system caused by the introduction of the correcting loops, 


Let us solve Equations (5) for the 1.a.c. of type I, The relationships needed for the other types of l.a.c. 
can be obtained from the corresponding Equations (5) in a similar manner, Therefore, the results for l.a.c. 
of types Il, Ill, and IV will be given below without proof, On the basis of (5) for a l.a.c., of type I, we have: 


os ~ BF (@)— Bra)’ 
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Because B(w) > 0, the minimum of oe at a fixed value of %c, 48 is apparent from (8), will occur when 
: B(w) is at a minimum, From Equation (9), after substitution 


, db . 
db deg of =28(y) - from (5), it follows 
W3 We D 
N that B(w) has a minimum at 
— 52) = (20) 
! Na According to (5) and (10) we obtain 
\ Thus, the minimum of “® occurs when the com- 
= Gh = pone 
la.c, of type I (see Equation (45) in the appendix) we 
\ ke % 
Fig, 2 Equations (19) to (12) establish for l,a.c, of type I 
the required relationship between the relative conjugated 
frequencies gy, and D, at the minimum value of the fre- 
quency interval ws. Hence: 
= — as 9) 


In the process of designing servosystems there may also be initially given one of the conjugated frequencies; 
but because the margin of phase stability in a system must in any case not be lower than 20°, it is not always 


e) possible to fulfill the above indicated minimum condition for the ratio ws/w, In this case Formulas (13) take 
on a somewhat different form, 
Let 
ae 
Then from Equation (5) for l.a,c, of type I we have 
or 
= (2— d)B(g). as) 


The relative conjugating frequencies are substantially positive quantities; therefore, Equations (14) and 
(15) make sense only provided the quantity d satisfies the inequality 


0<d<2. (16) 


) At d=1 the condition for the ratio w/w to become a minimum is fulfilled, According to Equations (12), 
(14), and (15), we finally arrive at: 
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Usually it is desirable to strive for fulfillment of the condition of having the smallest possible frequency 
band width in the l,a,c, with a slope of 6 db per octave near the crossover frequency at a prescribed margin 
of stability, that is to have d=1, However, ifd# 1, but has one of the values possible according to (16), it is 
easy to estimate how much the length of the asymptotic portion with a slope of 6 db per octave near the cross- 
over frequency is increased, 


On the basis of (13) and (17) we have: 


i 
(18) 


An analysis of Formula (18) shows that 1 changes slowly with d in the region of minimum 1 (d =1), 
For instance, for the interval 0.6 <d < 1,4, the increase in length of this portion of the asymptote in relation 
to its minimum length is not more than 19% (1 < {=< 1,19). However, when d is near zero or near two, the 
quantity begins to grow rapidly, At 0 <d <0,.3 and 1.7 <d< 2wehave = 2, 


Consequently, when designing servosystems, it is desirable that the quantity d should at least not be 
outside the range 0.4 < d < 1.6, 


The basic relationships used henceforth for determination of the parameters for correcting devices with 


1,a,c, of type I are given in Formulas (13) and (17) above, Corresponding expressions for l.a,c, of types Il, Il, 
and IV are presented in Table 2, 


TABLE 2 
Type of 1, a, c, 
Parameter 
I II Ill IV 
B (9) 8 (9) B (9) i 
B(¢) B (@) B (@) 
B (e) i B (9) 1 
Cc c 
i 2 2 
(2—d) B (9) | (2—d)8 () 


If the system includes links which have inertia and whose time constant Ty «< — (small parameters), 
c 
then the open loop transfer function of the System (2) can be represented as follows: 
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K(p) = K,(p) 


Il (1+ pT,) (19) 
ing 


Taking into consideration the smail parameters, as shown in the appendix, the equations that connect 
the relative conjugated frequencies and generalized parameters have exactly the same form as Equations (5) 


except that in the right side instead of 6 (g) there must be substituted g ,(y) , which is determined by the 
formula 


Ba = — Sr. (20) 
1—4 


Consequently the formulas of Table 2 can be used also for systems having a typical 1.a,c,, when small 
parameters are taken into account, if in the formulas of this table 8 (y) is replaced by Bn(¢), which depends 
not only upon the phase margin of stability but also upon the small parameters of the system, In addition, it 
is easy to show that the relationships that are already obtained, or that will be found in a similar manner, can 
also be utilized for servosystems with astaticism of the second order and can also be extended to systems with 
different open loop transfer functions [2]. 


It should be noted that the relationship obtained,as is evident from Table 2,are of a simple algebraic 
character, This makes it possible to calculate the parameters for correcting devices quickly, without laborious 
graphic computations, Besides, by means of Condition (16) it is possible to judge whether or not the designed 
system could be put into practical operation, 


2, Examples of Parameter Calculations for Correcting Devices 


1, Suppose the open loop transfer function of the given (uncorrected) part of the system is 


i 
it 
ima 


D i 


where T,,, is the time constant of the electromechanical components, T_ is the time constant of the ampli- 
fying and transforming devi::e, Dy is the amplification factor. 


Further let us assume that Tey, > Tx. It is necessary to choose the correcting arrangement in such a 
way that the over-all transfer function of the system shall correspond to a l,a,c, of type I and have the required 
values of we, and D, Thus, for a l,a,c, of type 1 (2): 


(22) 


it is expedient to use for the solution of this problem an integro- differential passive electrical network as shown 
in Fig, 3, and an inertialess amplifier. 
Indeed, assuming that 


D (1+ pT (4 + PTs) 


where D/D, is the amplification factor of the inertialess amplifier, Ty, Tg, Tem, and Tg are time constants of 
the correcting loop, Taking into account (3) and (21) we obtain 


D 1+ pT: i 
K, (0) = Kp) = (24) 
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After T,, is taken to be equal to Ts, and it is taken into account that Tg and T; are small parameters, 
a comparison of (24) with (22) shows that: 


K (p) = K; (p). (25) 


Thus, the problem is reduced to the determination of the unknown time constants T;, T3, and Ty, 
For the integro- differential network the following relationships 


hold true: 
(26) 
Al 
a— where 


= 2. 
Using the formulas of Table 2 we determine the time constants T, and T;, The sequence of operations 


n 
for this is as follows: from the given values of g and ®, >) 7, and the value for Tq just found, we obtain 
img 


4 
B, (9) = B(e) — > >) lem 
ima 


Then we determine is i and finally we obtain the t time constants 
hwe d=2 y sough 


After % and T; have been determined, it is possible to find the parameters of the electrical network from 
the following relationships: 
RC, =T—T,—Tem =Tem (28) 


For this purpose it is necessary to assume an arbitrary value for one of the sought for parameters, for in- 
stance the capacity C,, and then, using Formulas (28), find the rest of them. In assuming an arbitrary value 
for the first parameter, a standard size should be chosen, and moreover it should be chosen in such a way that 
the other parameters required to work with it should correspond or be near to standard sizes listed in a catalog. 


In the following a numerical example is given, and to illustrate conditions the open loop logarithmic 
characteristics of both amplitude and phase versus frequency are plotted * (Fig. 4). 


* All numerical examples and also the plotting of the logarithmic characteristics were carried out by L, E, 
Zalmanova and A, A, Vlasova, 


om | | 
But because (see appendix) we have: 


Given =12sec *, D = 300 sec™*, Tey, =0.1 sec, Ty = 0,05 sec, Ty =0,01 sec and = 30°, 


We take T,, to be equal to Ty = 0,05 sec (ws = 20 sec”); then Bn(g)~w 0.44, d= 0,64, which satisfies in- 
equality (16), T, 0.27 sec ( = 3.76 sec™*); ¥6.67 sec (uy = 0.15 sec™*), 


de Now let us determine the parai.seters of the correcting 
loop, Let Cy = 10 uf; then Ry = Tz/C, = 27 kohm, 


RyC, 6.3 sec, Ry = 630 kohm, 


" The resistor nearest to 630 kohm in the catalog is 620 kohm, 
Substituting this value of Rg into RgCy = Tem = 0.1 sec, we obtain 
Cy w 0.16 wf, We finally set Cg=0.15 pf, 


30 

: sod A 2. Suppose that the transfer function for Ky(p) is given in 
wh the form of (21), and the time constants T,,,, and T;, are close 


log to each other and both satisfy the inequalities T,,, < ram and 
c 


< It is required to choose correcting device in such 
20 “i \ a way that the transfer function of the system shall correspond to 
T 
Fig 


a l.a.c, of type Ill, and that we, @c, and D attain given values, 
On the basis of [2], for a l.a.c, of type Il we have 


= 1+ (29) 
K 


then it is easy to find, that 
Di+ 
that is, the correcting device represents an inertialess amplifier with an amplification factor of = and a 


passive integrating loop, as shown on Fig, 5. K(p) will equal Kyyy(p) if Ts is made equal to V¥Tem Tk and 
Tj is taken into account as a small parameter, 


k, The problem then is reduced to the determination of T; 
o———{___}—_++ -~ <2 and Ty, The sequence of the procedure is the same as in the 
first example, except that d is found by solving the equation 


% for = given in Table 2, Column Ill, which yields: 


d=2 | 
4 
Fig. 5 The parameters of the integrating loop are determined 
by means of the formulas: 


T, = RC. 


Here is a numerical example. The following data are given: = 5 sec™*, D= 200 sec™!, = 30°, 


Tem 70.08 sec, Ty =0,07 sec, Tj =0,01 sec, We calculate Ty to be equal 0,075 sec of wy *13,3 sec™* and 
then find: 


n 
B, (*)and (9) = B(9)— = 0.5, d= 0.5. 
* 
Then T, and T, and also the parameters of the integrating loop are determined: 


Bm 
ar 
| 
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j 


=0.78sec = 1.28 sec" ), = 31.3 sec 0.032 sec 
c 2 


Let C=20 yf, then Ry = 39kohm,R, =n (7-1 ) 1,6 megohm. 
On Fig, 6 are given log log plots of amplitudes and phases versus frequency for the open loop system of 


this example, 
d 
4 
60 
2 
\ 
/ 
\ yi 
y 
Ni 
20 \ 
Fig. 6 
% 
Fig. 7 


3, Let K,(p) have the form of (21) and let 1/Tem < 
< We. It is sought to choose a parallel correcting device so 
that the over-all logarithmic amplitude characteristic of the 
system will be of type Il at given values of we, ¢¢, and D, 


We throw around the link that has the highest time con- 
stant Ter a flexible tachometric (velocity responsive) feed- 
back line and include in this line a two section high-pass RC— 
type wave filter (Fig. 7). 

Then we have: 

Koo (P) = 
Koop 
1 
[1 +20; + (31) 


(where Tr is the time constant of the filter), or when ~ « 2, 


ph + (32) 


(P) = 


Comparison of Equations (4), (21), and (32) shows that 
the over-all l,a.c, of the system is reduced to type Il if we put 


Ts =T2, Ty =Ts, 
(A+ pT, + + = (1 + PTs) (1+ PT (33) 


where T, is a small parameter. 


Equating in (33) the factors by which p and p” are multi- 
plied on both sides of the equation we obtain two equations, 
from which it follows that: 


Rkmp , 
D 2 D 
V2 r,). (34) 


Thus, the problem now is reduced to finding T, and K,,. 
The procedure for this is as follows: using the initially pre- 
scribed values for we, @, and Tj and the value for Tg cal- 


culated from (34), we find 6 ,(y) by means of the formula: . 


img 


Then we determine the value of d from d=2 — me 
numerical values of 6,,(g) and d, we find: 


» Where ws = = , and finally, knowing the 


@) 


and then K,, by means of Formula (34), 
The parameters of the RC filter are determined by means of the relationships: Ts = Tg = RgCg = RyCy. 


We start out with taking as given one of the parameters and the relationship ag « 2 and then de- 


termine the rest of the parameters, 
We now give a numerical example, The following 
data are given: we = 20 sec™’, D = 400 sec™*, = 
=0,1 sec, =0,025 sec, Ty = 0,002 sec and = 30°. 
spre. We find Ty = 0,005 sec, d= 0,89,8n()~ 0.45, 
0.195 sec, Ko, = 18.2, 


Then for the parameters of the RC filter we have 


\ = RyC, = 0,195 sec, Let Cy = Quf and Ry/R, = 0,1, Then 
Ry = 97 kohm, Wetake Ry = 100 kohm, and we then obtain 


Rg = 1 meg. and Cg=0,2 wf, 
t On Fig, 8 are shown the open loop logarithmic am- 


T ‘ < 7 plitude and phase versus frequency characteristics for this 

example, 

| | 4. We are given the open loop transfer function of 
a 
1 Tla+pr,) 

ot where T;, and T, are time constants close to each other, 


1 1 
Fig. 8. both smaller than ——,and T,,..> 
ig We om We 


It is required to select a correcting arrangement in such a way, that the over-all l.a,c, at the given 
values Of we, Mc, and D shall be of type IV. 


C In solving this problem, we can utilize an inertialess 
| amplifier, a parallel correcting device, which has a transfer 
Q function of the form of (32) and a series correcting device 
R, which has a transfer function: 
R, 
k (P) (36) 
where y is chosen from the condition Tyy < Ts, This 
- ‘ “= transfer function can be obtained by means of a differentiat- 
Fig. 9 ing electric network (without taking into account the loss) 


as shown on Fig, 9. The system is then brought to type IV 
of the l.a.c,'s if we make: 


T + P+ + Kee PM = (1 + (1 + pT,) 
and account for the time constants Tg and T,y by considering them as small parameters, Thus, the problem is 


reduced to finding of T, and Kg, and to the choice of y. The values of T, and K,, are determined in exactly 
the same way as in the preceding example; only, the value of d must be found by means of the formula: 


| | 
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@=2(1- 


The amplification factor of the inertialess amplifier is chosen to equal Da 
Let us now introduce a numerical example: Given 
We =18 sec™*, D= 400 sec*, = 30°, Tern = 0.1 sec, 
Ty = 0.02 sec, T,, =0,.015 sec, Ty =0,03 sec, Tj = 
= 0,001 sec, 


We find Tq =0.0045 sec, Ty = 
= 0,017 sec (ws = Let y = 0,2; then 


ine By = B (e) — 
-+ (Ty + Ta +Tyy) = 0.42, 


aX 

d=0,53, T;= @, (9) 

=0.395see, = 2 —1— 


| 

| The parameters of the two-section wave filter are 
ag 4 selected in a way similar to the one used in example 3. 
Let us find the parameters of the differentiating network 
\ shown on Fig, 9. RgC = T,; we set C =0.05 wf, then 

Ry =600 kohm; we take Ry = 620 kohm, Rg= y RyAl— y) 
is calculated to be 155 kohm, We take Rg = 150 kohm, 


Fig. 10 The corresponding logarithmic characteristics are 
plotted on Fig, 10, 


In a similar manner with respect to the basic idea as well as to details of procedure, the parameters of 
correcting arrangements can be determined for other uncorrected K,,(p) functions, 


<= 8 
L 


APPENDIX 
Proof of Equations (5) ’ 
It is known that: 
% =F + (37) 
K 4, = 1, (38) 


where K(w) and gw) are respectively the amplitude versus frequency and the phase versus frequency character- 
istics of an open loop servosystem. 


For l,a.c. of type I according to (22), we have: 


9 (oe) =— —arc tan tare —arc tan (39) 


For the frequency band wz <w <ws (in the neighborhood of the asymptotic crossover frequency w,), the 
following inequalities hold true: 


| 
| 
Ix 
— 
tt 
| 
| 
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a a 
We expand all three arctangent functions arctan ~~ (39) into series, and then use the approximation 
obtained by discarding all members higher than of the first order, Equation (39) then is changed into: 
(41) 


After substituting (41) at w=w, into (37) and rearranging, we obtain: 


S (1-3) 20, (42) 


The ratio rs as a function of we, and D is obtained from Equations (22) and (38), From (22) 


iter (43) 
Taking into account (40), Formula (43) can be simplified into 
K = 2 = (44) 
According to (38) and (44) we have: 
 % 
(45) 


Finally, after substitution of (45) into (42) the latter becomes identical with (5): 
(4 ba 


In a manner similar to the above, equations corresponding to (5) can be obtained for 1,a,c,"s of types I, 
Il, and IV, 


Accounting For Small Parameters 


For transfer Function (19) the phase versus frequency characteristic is: 


(@) = (@) — Djare tan”, (46) 
ima 


where g(w) is determined by Equation (41), 


Because T; = +. <- , the inequality « 1 holds true and we have: 
Wi 
n n n 
ra) Mo 
img 
We then substitute Formula (46) putting w = we into Equation (37) and obtain: 


n 
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Introducing: 


we obtain an equation of the same type as (5) except that on the right side  (y) must be replaced by 6p(¢). 
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AN ANALYTICAL FORMULATION OF THE SYNTHESIS PROBLEM OF 
CORRECTIVE DEVICES IN LINEAR SERVOSYSTEMS 


Vv. G, Segalin 


(Moscow) 


An analytical formulation of the generalized problem of synthesizing corrective 
devices in linear servosystems is proposed, and, on the basis of the initial synthesis 
equations which are derived, a method is worked out for determining the transfer func- 
tions of corrective devices for the most frequently encountered cases of correction, In 
conclusion, an example is given of a solution of the synthesis problem, 


INTRODUCTION 


At the present time the least developed branch of the theory of linear electrical circuits with lumped 
constants and linear automatic control systems is the theory of their synthesis, including the most important of 
its branches — the theory of synthesizing corrective devices, 


Notwithstanding a series of major achievements in this branch which have been reported in [1-6] and in 
a series of other works, there is lacking up to the present time: 


1) An analytical formulation of the synthesis problem for corrective devices; 


2) A method of determining the optimal transfer function, which the corrective system's transfer function 
should be made to equal, and its generalized analytic expression; 


3) A formulation of recommendations for the selection of various corrective networks, keeping in mind 
the peculiarities of the invariant parts of the servosystem and the actual uses to which it is put; 


4) A formulation of recommendations as to the expediency of using active corrective devices, again 
keeping in mind the actual conditions under which they will be used, 


In the present paper an attempt is made to formulate the general synthesis problem analytically and 
generalized transfer function equations for corrective devices are derived subject to the conditions that their 
expressions take into account requirements of type [7] for the corrective servosystem as well as the peculiarities 
of its invariant part, 


1, Formulation of the General Synthesis Problem of Corrective Devices 


Let us take some open linear electric circuit, or open servosystem, or automatic control systerm made up 
of N directed sections, where the transfer function for the i-th section is known to be 


(Bo, By, By + BaP") 


Wi(p) = 


The entire open circuit's transfer function is 


N 
Co+C Cop? +... 


: (2) 
im} p" (Do+ Dip + Dep* +... + 


Let us assume that this servosystem is not functioning properly and that m special corrective quadripoles 
(CQP) are inserted into the corrective circuit (C-circuit) to correct this, 


It is also assumed that the transfer function (2) corresponds to the invariant part of the servosystem and 
that the only variable parameter in this servosystem is the resulting constant gain factor, Consequently, it is 
also possible to vary the constant gain factor of separate sections,* 


On the basis of the foregoing, the transfer function of the C-circuit after correction is 


(Co + ACo) + + p+ ...+(C,, + AC,,) p™ 


ph + + (Di + p+... (Dy + AD,) (3) 
N 
I] Wi(p) Up 


which we will call the desired transfer function (GF) because it is the result of the corrections, 
In Expression (3), we designate 


Keg= (Ko + (et 


as the new constant gain factor of the servosystem, and + AKy, +ACy, +ADj, + AF, +AC; as the increments 
added to the corresponding constant coefficients of the C-circuit's transfer function as a result of the correction, 
while S (p) and U(p) are polynomials in p: 


S(p)= 1+ 
U(p) =1+U p+... Uypr. 


It is assumed that since the synthesis problem has been solved, it is possible to determine the set of 
COQPs, introduced into the C-circuit, the plan for their insertion, and the analytic expression for the transfer 
function of each of the CQPs, 


It is assumed further that the corrective networks selected, and the set of CQPs inserted, make it possible 
to get the simplest analytic expressions for their transfer functions, 


It is assumed further that: 


1) The introduction into the circuit of active corrective devices (CQPA) is undesirable because of their 
relative complexity as compared with passive corrective devices (CQPP) and can be justified only in those cases 
where the use of CQPP does not yield the desired results; 


2) Since the method of determining the circuit of a CQP, given its transfer function, is known [8-10] 
and will not be considered in the present paper, the solution of the synthesis problem in the present paper will 
lead to the determination of the corrective device's transfer function, As far as imposing realizability re- 
quirements on the expressions we have calculated for these transfer functions is concerned, this matter will be 
dealt with in another article, 


* Although this limitation complicates the proposed problem, it also makes it possible to solve the synthesis 
problem for such systems, the invariant part of which consists of standard assemblies and components, where in 
most cases the only variable parameter is the constant gain factor, 


2, Networks For Introducing Corrective Devices 


To get exhaustive recommendations as to the expediency of using various corrective networks, it is 
necessary, obviously, to consider the greatest possible number of variants of these networks, 


In the author's opinion, the networks for introducing CQPs into C-circuits which are listed below and are 
depicted in general form in Figs, 1-9, include the overwhelming majority of corrective networks known to the 
literature, As far as the nomenclature for all the varieties of corrective networks that are listed is concerned, 
it is quite arbitrary and is introduced solely for convenience of presentation, 


= 
Synth | $ 


Mm(p) | ~ Oyom(P) 
=. | { 
M,(P) 
iD, 
M(P) Wf) > 
ir j 
4 } 
ases 
} 
ll 
me Fig, 1 Fig, 2 Fig. 3 


Appendix I gives solutions of systems of operator equations which correspond to all the corrective net- 

; works depicted in Figs, 1-9, The analytic expression for Ky (p) - Kg(p) which is calculated there is the transfer 
in function for the C-circuit, modified by the insertion of the indicated network of corrective devices, In these 
expressions the unknown transfer functions for the CQPs which have been introduced appear as the operator 
functions M;(p) which we try to determine, where i (i=1,... , m) corresponds to the number of the corrective 
device, 
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3, The Generalized Equations of Synthesis of Corrective Devices 
and Their Solution 


In accord with the above the most general form of the synthesis equation is 
Ko(p) = Kg (p), 4) 


where K,, (p) is the transfer function of the corrective circuit after the insertion of m corrective devices 
according to any of the networks depicted in Figs. 1-9 (see the corresponding expression for Ky(p)+Kg(p) 
in Appendix I), o is the number (g =1, . . « , 9) of the corrective network corresponding to Figs. 1-9, Kg(p) 
is the desired transfer function, in general form corresponding to Expression (3), 


} 
(P) 
a + 
lng?) 
Wysm(P) Mm(p) 
| 
! 
A 4 4 A 
| 
Wag > ~ 
> 
i 
My(p) W, (2) M,(2) 
> = = 
| 
W,(p) | 
> > 


Fig. 4 Fig. 5 


We will recall that in Equation (4) and its particular cases, the transfer functions Mj(p) are the ex- 
pressions sought, 


Particular cases of the synthesis Equation (4) are considered in detail in Appendix I, corresponding to all 
nine types of corrective networks (Equations I-IX) considered in the present paper, 


Pout? 
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The solution in general form of the above-mentioned equations makes it-possible to determine the 
analytic expressions for the corrective devices’ transfer functions which in this case will have the following 
form, 


A, For series correction networks 


il M (p) S(p) 
(5) 
Ww, (p) 
Ws 
Oy. m(?) 
AE 
(P) Mm 
{ges 
z 
Wy. M,(?) 
| 
M,(p) 
he 
Wy, 
Fig, 6 
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B, For parallel correction networks 


I, Cases of the directly-connected type: 
1) Method of bypassing one at a time 


Y,(p) 
Mi (p) = (i= 1,...,m); (6) 


2) First type of step-by-step correction 


i 
M; — ™ 
i (p) Weis (p) m) (7) 


3) Second type of step-by-step correction 


Y; (p) 


a;Z(p) Il (i= 1,...,m); 


jei 


M; (p) = 


Wy, 


Wary (p) My) M,(p) 


> 
AA 


Fig. 7 
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4) Method of cross-switching 


(Pp) 
j=1 
(0) 
(2) 
| Mmm (0) 
> 
4% 
| 
| M,(p) 
By 
| 
(P) M,(p) wif 
W,(p) 
Fig, 8 
Il, Cases of the feedback type: 
1) Method of bypassing one at a time 
¥,(P) 
2) First type of step-by-step correction 
Y;(P) 
== 
a,2;(p) ™ (11) 
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3) Second type of step-by-step correction 


Y; (p) 
Mi (P) = (i= 4,...,m); (12) 
II (Pp) 
4) Method of cross-switching 
Y; (p) 4 
M1 (Pp) = m (i= 1,...,m). 
i4; (P) (13) 
II Ww; (p) II (Pp) 
jimi 
It follows from consideration of the equations we have 
f derived that the solution of the synthesis problem in the 
Wom(0) form we have proposed becomes the determination of cer 
T tain intermediate operator expressions 
= Yi(p), Zi(p) 
Woms(D) My, (2) 
Zi(p)’ (p)’ 2; (p) 
— al It is possible to calculate all the unknown constant 
: * coefficients of the polynomials Y;(p) and Z; (p) from 
equations which are given in Appendix I, 
Mns(P) 
> | 4, Discussion of the Results 
= 


ae The fundamental result of the present paper is the 
' derivation in general form of analytic expressions for the 
transfer functions of corrective devices, 


- — __ | Practically all corrective device networks of the series 
Wm (PD) M,(p) connected, direct connected, or feedback types known to the 
literature were considered, The methods developed also may 
Pe + be extended without difficulty to such types of insertion of 
corrective devices as were not considered here, 
Won (p) 
The following limitations are imposed on the class 
~ of servosystems considered: 
Wrre(P) 1, The corrective servosystem and the corrective 
Pa devices inserted into it should be made up of linear sections 
<= it 
with lumped constants, 
' } 2, The sections inserted into the servosystem should 
W,(p) | | be detecting (directed action), 


Rane As we see, the limitations listed are minimal, They 
- do not prevent consideration of a wide class of the most 
3 frequently encountered linear automatic control systems and 
servosystems, made up of circuits with lumped constants, 
Fig. 9 and linear electrical circuits in which wave phenomena may 
be neglected, assuming, that they were made up entirely of 
lumped constants, 


Passing over to a discussion of the problem of actually constructing equipment corresponding to the 
corrective devices" transfer functions, it is necessary to note the following, 
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All the expressions we have gotten for the transfer functions discussed may be constructed as active or 
passive CQPs or combinations of both, depending on the correction problem involved, 


Analyzing the expressions we have gotten, we may also note that imposing the requirements listed on 
these transfer functions in its turn makes necessary certain restrictions on the structure of the transfer functions 
which make corrections in various sections of the C-circuit, In this way, requiring passive CQPs or imposing 
other conditions on the corrective devices is not always compatible with the peculiarities of the invariant part 
of the servosystem, On the other hand, the expressions found for the CQP's transfer functions, when these 
corrective devices are made up of passive elements, have mutually contradictory structures when these corrective 
devices are inserted as a directly-connected or feedback network, We have in mind the differences involved in 
satisfying the requirement that the degree of the polynomial in the denominator be greater or equal to the 
degree of the polynomial in the numerator of these transfer functions, From this, it is possible to make a far- 
reaching conclusion, amounting to the following, 


By virtue of the mutual incompatibility of the limitations connected with the structural peculiarities of 
the transfer functions of sections of the corrective servosystem, the correction methods by means of directly- 
connected and feedback networks are mutually supplementary in the sense that the impossibility of using one 
method because of the conditions imposed becomes the optimal situation for applying the other method, 


Upon analysis of the expressions we have derived we come to another no less important conclusion, 


The sign of the CQP's transfer function is determined by the sign of the coefficients in the intermediate 
expressions 


Y; (Pp) Y; (Pp) 
Yi(p), Zi(P), 1— Z; (P) Z,(P) — letc, 


that is, the sign is determined ultimately by the signs of the increments, added to the time constants in the in- 
variant part of the servosystem as a result of the correction, Therefore, when the question is put in general 
terms, it is necessary to remember that in using parallel corrective networks,* each insertion of a CQP into 
the corrective circuit can have either a positive or negative sign and that this is determined by the sign of the 
analytic expression of the CQP's transfer function, Therefore, we must mention that, in general, CQP"s may 
be inserted in a C-circuit by means of positive or negative directly-connected networks, or by positive or 
negative feedback networks, 


Analysis of the expressions we have determined also leads to the following generalizations, 


1, Since the constant gain factors of the various sections have a direct relation to the absolute magnitude 
of the coefficients of the polynomials Y;(p) and Z;(p), and consequently, on the signs of their differences in the 
expressions 

Y; (p) Y; (p) 


etc,,then by suitable changes of the constant gain factor of various sections of the invariant part of the servo- 
system** it is possible to make the desired changes in sign of the connection when different corrective devices 
are inserted (for example, substituting positive feedback for negative), 


2, Only through a suitable and therefore unique choice as to how to connect the corrective network 
(directly-connected or feedback) does the proposed synthesis problem become solvable when it is required that 
the analytic expressions of the corrective devices’ transfer functions, given a certain invariant part of the 
servosystem and a certain desired transfer function, correspond to the use of passive corrective devices, However, 
the use of active corrective devices in solving a practical synthesis problem may permit the partial, and in some 
cases, the complete relaxation of the structural limitations which were established here, and which must be 
applied to the transfer functions of sections of the servosystem's invariant part or to the corrective devices’ 
transfer functions, Nevertheless, because of the mutual incompatibility of the limitations associated with the 


* In one given servosystem, 
** See footnote on p, 138, 


structural peculiarities of sections of the servosystem's invariant part, when using directly-connected or 
feedback type networks, use of CQPA is usually inexpedient and can be justified only in those rare cases where 
the corrective network is subject to severe restrictions and the correction problem can not be solved using only 
COQPP, 


In the present paper, the problem proposed was to find general analytic expressions for CQP's transfer 
functions for the broadest class of networks that are inserted into C-circuits, subject to the condition that GF 
also be given in general form, 


The next step in the investigation is the establishment of all those requirements which must be satisfied 
by the servosystem's invariant part and by the CQP's transfer functions in order that the synthesis problem can 
be solved in the simplest manner, 


APPENDIX I 


1, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig, 1, 
The initial system of equations is 


6: (Pp) = (p) (P), 
8; (p) = (p) (p), 
(p) = 8; (p) Ws(p), 


6, (Pp) = % (p) Ms (p), 


(P) = Om (P) Wn (P)s 


out (P) = On (P) = (P) Mmm (P)- 


The unknown transfer function is the solution of the following system of equations: 


im] 


K;(p) = 


The synthesis equation is derived by substituting in (4) the expressions K,(p) and Kg (p) and getting 


WwW M,(p) = 


The solution of this equation with respect to the m unknown corrective devices’ transfer functions gives 
Expression (5) (see text of article). 


2, Find the C-circuit’s transfer function when corrective devices are inserted as in the network in Fig. 2, 
The initial system of equations is:* 


(p) = 6, (p) (P), 
(P) = 210, (p) Mi (p) 


(Pp) (p) + (P) (Pp) aM, (p)| 6, (Pp), 


(P) = 04 m—1 (P) Wy (P)s 
* Here, and below, with the exception of the cross-switching cases for the CQPs, it is assumed that the corrective 


devices are inserted in the C-circuit somewhere between the v+1 and v+m _ sections, and that the coefficient 
a; gives the fraction of the signal that is found in the given part of the C-circuit, 
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8, 4m (P) = (P) + (P) = [Wy 4m (P) + My, (P)) Oy 


N 
Gout (P)=O4m(P) Wi 


The unknown transfer function is the solution of this system of equations 


N m 
a,M, (p) 
K. = | w | | 
2(p) i (p) Ww (p) 


By substituting K,(p) and Kg (p) in (4) we get the synthesis equation as 


N 


i—1 


Expression (6) for the i-th corrective section may be found by solving the above equation, To find the 
unknown coefficients of the polynomials Y;(p), 2Z;(p) in Expression (6) we consider the equation 


which in expanded form is 


fl = I 3, (p) S(p), 


i=1 i=1 
where we assume 


4, (p) (p) 
Wi P= MiP) = 


4, (p), (p). 3; (p). 3; (p) are polynomials in p. 


3, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig. 3, 


The initial system of equations is 
6,4, (P) = 0, (p) (p), 
0,41 (P) = 410, (p) Mi (p), 
8,41 (P) = (P) + (Pp) = IW, 4 (p) + 8, (P), 


|| 
ve 44m (Pp) (Pp) (Pp), 
t 
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(P) = 0, (P) M yn (P), 
(P) = (P) + (P) = Wy (P) 9, (P) + mm (P) 9, (P), 


0, (p)= 0, 


i=1 


Gout (P)=%4m(P) 
i=—v+m-+1 


The unknown transfer function is found by solving the following system of equations: 


N m 
M 


By analogy with the foregoing the synthesis equation is 


(Pp) S (p) 
j=1 


The solution of this equation gives Expression (7) for the i-th corrective section, The unknown coefficients 
of its polynomials Y;(p), Zj(p) may be found from the equation 


¥; (P) Wi(P) s(p) 
(P) 


which may be written in expanded form as 


m i-1 
DM Il Z;(p) = Il 4; (P) 4, (P) I 3; U (P), 

Tl 3; (p) U (p). 

i—1 i=—1 


4, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig, 4. 
The initial system of equations is 


(p), 


i= 
(p)= (p) a, (p)6 v+m—1 (p) = 


= Sam, (p) 6,4 


8, 44-1 (Pp) = (p) 


i-1 


jel 
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0,(p) = (p) 
N 


|] Wyle). 
i=mv-+m+1 


The unknown transfer function is found by solving the following system of equations: 
N m a,M, (p) 


By analogy with the foregoing the synthesis equation is 


N 
i~1 — Wes; (p) i=1 


The solution of this equation gives Expression (8) for the i-th corrective section, The equations for finding 
the unknown coefficients of the polynomials Y;(p), Z;(p) are analogous in general form to the previous case, 


5, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig, 5, 


The initial system of equations is 


(p) = 0, ,(P) (P). 
i=] 


(p) = 0, (Pp) Mi (), 
(Pp) = (P) + (P), 


(P) = —g (P) W om —y (P)s 
—1 (P) = Om (P) (P), 


The unknown transfer function is found by solving the following system of equations: 


i—1 
im1 


2m=N m 
T] Wp) | 1+ 
- 


By analogy with the foregoing the synthesis equation is 


m 2m=N 
- 


Il W; (Pp) (Pp) 
j=i 


(v) 


The solution of this equation gives Expression (9) for the i-th corrective section, The equations for finding 
the unknown coefficients of the polynomials Yj(p), Zj(p) are analogous in general form to the previous two 


cases, 


(p) 
| 
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6. Find the C-circuit’s transfer function when corrective devices are inserted as in the network in Fig. 6, 
The initial system of equations is 


0,41 (P) = 10, (p) + (Pp), 4, Wy, (P), 


0, (p) = 0, Wi 
i=1 


out (P)=%4m(P) [] Wy (P). 
i—v+m+1 


The unknown transfer function is found by solving the following system of equations: 


K = 


N 


Ww S (p) 
Tl Il Wi () U(p) (v1) 


The solution of this equation gives Expression (10) for the i-th corrective section, To determine the un- 
known coefficients of the polynomials ¥;(p), Z;(p) in (10) we consider the equation 


m N 
Z; (P) WilP) S(p) 


which may be written in expanded form as 
m N N 
im1 


m N N 
I] ¥: =T] 3; 


i=1 
7, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig. 7. 
The initial system of equations is 


(p) = 6, (p) + 21M; (p) (P), 


0, (p) = %in(p) |] 
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6,44 (P) = 10,44 + a,M i (p) (p), 
(P) = 10, (P) + My (P) 9, (PI (PD 
N ; 


N 
our (P) = (p) Il W, (p). 


The unknown transfer function is found by solving the following system of equations: 


i 


N 
Ky (p)-= |] 
ft 1 — aM, [] W447) 


im1 j=i 
By analogy with the foregoing the synthesis equation is 


S(p) 
Il (Pp) m m = Ik Wi (Pp) U (p) (VID 


The solution of this equation gives Expression (11) for the i-th corrective section, To determine the 
unknown coefficients ef the polynomials ¥;(p), Z;(p) in (11) we consider the equation 


¥ N W,(P) U (p) 
27, (P) Wi (p) 


which may be written in expanded form as 


m i—1 m 
DU 
i—1j=—1 jmi+1 
N ‘ N N N 
i—1 ie 
m N N 
IL 
i—1 i—1 t—1 


8, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig. 8, 
The initial system of equations is 


(p) = 0, (p) |] War), 


8, (p) = 8 (p) + 05° (p), 
6 

Il 


[] 
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j=i | 
v 
m 


N 
Gout (P)= %4m(P) Wale). 
t=—v+m-+1 


The unknown transfer function is found by solving the following system of equations: 


K,(p) = [] Wi (2) i 
1— aM, (p) T] 
im] j=) 


By analogy with the foregoing the synthesis equation is 


S (p) 
I] W, (P) i 


j=1 


(vi) 


The solution of this equation gives Expression (12) for the i-th corrective section, The equations for finding 
the unknown coefficients of the polynomials Y;(p), Z;(p) are analogous in general form to the previous case, 


9, Find the C-circuit's transfer function when corrective devices are inserted as in the network in Fig, 9, 


The initial system of equations is 


8in(P) = %,, (P) Ma (P), 
(P) = + (P). 


(P) = (P) Mp, 

0,,, (Pp) = 95,_, (P) Wy (P), 


(P) = out (P) = %m(P) [] Wile) 
i=m+1 


The unknown transfer function is found by solving the following system of equations: 


2m=N 


i—1 j=1 j=1 


By analogy with the foregoing the synthesis equation is 


4 S( 
Il VW, (Pp) ™m Il Wie) 


i— > a,M, (p) (Pp) Il Wj (Pp) 
imi j=t 


(IX) 


The solution of this equation gives Expression (13) for the i-th corrective section, The equations for finding 
the unknown coefficients of the polynomials Y;(p), Z;j(p) are analogous in general form to the previous two 


cases, 
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APPENDIX Il 


An example of a solution of the synthesis problem 
Given a servosystem with first-order instability whose transfer function is 


| 
K (p) Ws (p) Ws (Pp) Wap) We (p) Wo), 


where 
Ky = 100 
x, + Uap) _ +6.25p) (1 + 0.6679) 
Walp) Ky + Uap) _ g(t +0.74p) (1 0.008679) 
2(P)=Ks (1+ Tap) ~ +0.952p) + 0.005p) 


i i 


i 1 i i 
WalP)= Tp) = Tp) = TF 0.0008) 


The correction problem is reduced to the following. 


It is necessary to find: numerical formulas for the transfer functions of the corrective devices which will 
be used in feedback networks; their number; where they will be inserted in the servosystem, etc, As 4 result 
of the correction, the servosystem should satisfy the following requirements: 1) have first order instability; 2) 
after the completion of one step of step-by-step control it should not have an overcorrection g. exceeding 25%; 
3) it should not have a transient time exceeding 0.85 seconds after the completion of one step of step-by-step 
control, Finally, the constant gain factor of the entire open servosystem and its various branches should remain 
unchanged after the correction, 


Using the methods described in [1], it is possible to calculate the optimal transfer function, which the 
servosystem should have after correction: 


Ko (1 + Usp) (1 + Usp) (4 + Usp) (1 + Uap) 
100 (1 4+ 6.25p) (1 + 0.667p) 
+Tep)d (i + 0.003p) (1 + p) 


(1 +0,74p) (4 + 0,00667p) 
X 0.004p) (1 + 1.56 p) (1 + 0.00125p) (1 + 0.0009p) * 


where AT, = 2,2 seconds, AT, = — 0,532 seconds, AT; = 0,048 seconds, AT, = — 0,001 seconds, 


Turning to the choice of corrective device circuit diagrams, we note that the following conditions are 
implied by Expressions (10) - (13), In order to be certain of an exact solution of the synthesis problem while 
using passive corrective devices (not modifying the optima] transfer function by inserting additional inertial 
sections), it is necessary to include in the correction scheme only those sections whose transfer functions have a 
numerator whose order is not less than the denominator's order, 


To satisfy these correction requirements only the first and second sections can be used, In order to get 
the simplest network of corrective devices, we will insert them into the corrective servosystem by the method 
of bypassing one at a time,* 


From the above, and using (10) and (VI), we get 


* A detailed comparison of the correction circuits considered in this paper will be made in another article, 
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(11 +475) ph A72) + (7s + 47s) +1 Ph _ 
(1 + Tip) (1 + Tap) (1 + Tsp) (1 + Top) 


(1 + 10p) (1 + 0.003p) (1 + p) (1 + 0.004) 
(1 + 7.8p) (1 + 0.535p) (1 + 0.952p) (1 + 0.005p) 


and further 
Zs (p) = 1+ 2, p+ 2, op? = (1 + 7.8p) (1 + 0.535p), 


Z2(p) = 1+ Z, oP + = (1 + 0.952p) (1 + 0.005p). 


In contrast to the way the polynomials Z; (p) were chosen, it is necessary to define the polynomials Y; (p) 
so that the difference Z;(p)— Y;(p) should have coefficients all of the same sign, Only one case satisfies this 
requirement: 

Yi(p) + = (1 + 10p) + P), 
Y2(p)=1+Y,,P + = (1 + 0.0033p) (1 + 0.004p), 


since in this case the difference of the polynomials 


(p) — ¥1 (Pp) = — P + (21,9 — ¥y,9) P? = — (2-66 p + 5.83p*), 
(Pp) Ys x 21) P + 21) = 0.95p + 0.00475 p? 


does not have coefficients of different sign, 


Using this method to calculate the unknown coefficients, and Formula (10), which is relevant here, we 
compute the unknown formulas for the corrective devices’ transfer functions: 


Yi(p) 


p(1-+2.19p) 
= Weil! — Zp) 


(1 + 6.25p) (1+ 0.667p) ° 
1 Ys (p) p (1 + 0.005p) 
Since the method of proceeding from the expressions we have derived to actual electrical networks has 
been worked out in sufficient detail, we consider that the solution of the proposed synthesis problem has been 
completed, 


We confirm the correctness of the result by substituting the expressions we have computed into the formula 
for Kg (pj), Since Kg(p) is found to equal the optima! transfer function, we conclude that the synthesis problem 
was solved correctly, 
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In conclusion, we note that having computed the corrective devices’ transfer functions, the signs of these 
functions indicate that the first should be inserted into the servosystem as a negative feedback network, the 
second as a positive feedback network, 
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AN AC BRIDGE WITH SELF-BALANCING IN TWO PARAMETERS 


V. iu. Kneller 


(Moscow) 


An ac bridge with self-balancing in two parameters in the same arm is considered, 
The effects of interconnection between the self- balancing circuits and of nonlinearity in 
the static characteristics of the bridge on the stability are analyzed, A bridge circuit in 
which it is possible to prevent interaction between the self-balancing circuits over a 
wide range of measured quantities is presented, 


1, Introduction 


Two time-varying components of an impedance (R and X, C and tan6, etc.) have often to be measured 
accurately in automatic process controls, 


An ac bridge with self-balancing-in two parameters is the best device for the purpose, The circuits for 
this so far given can be divided into two classes, according to the method of producing the signals which control 
the two driven elements, 


__ In the first group, described in [1-3], the control signals are produced by splitting the output voltage from 
the bridge diagonal into two components using two phase-sensitive null detectors, Figure 1 shows the block 
diagram of one such bridge, In general both control signals (€, and €3) are functions of both controlled para- 
meters, and hence the balancing circuits are coupied, 


~ 


4 


Fig. 1, Block diagram of a bridge, M) bridge circuit; P) preamplifier; 
N) phase-sensitive null detectors; SA) servoamplifiers; SM) servo- 
motors; Uy, and U,) reference voltages of phase angles y and @ re- 
lative to I. 
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In the second group system proposed by Karandeev et al, (3, 4] the control signals are produced by using 
two special differential indicators (differential balancing detectors), Each detector is supplied via phase-shifting 
networks with two voltages taken from different points in the bridge, The detectors produce signals which mainly 
depend on one parameter only, which means that it is possible to prevent interaction between the self- balancing 
circuits over a wide range of measured quantities, 


However, these detectors cause a considerable increase in the static error (up to 2%) (5), Karandeev and 
Grinevich [6] have effected some reductions in this source of error, None the less, it is very much more difficult 
to design these detectors for stability and high precision than in the case of the first group of circuits, 


The first- group bridges have only been analyzed [25] in relation to the measurement side, It has there- 
fore been (incorrectly) supposed [3-5] that the coupling between the balancing circuits and the nonlinearity of 
the static characteristics will mean that the first group will in principle have unsatisfactory transient response, 


The above argument shows that the first group is of interest, and that an analysis of the transient response 
is ot value, The method of producing the control signals using the phase-sensitive detectors is equally applicable 
to complex automatic balancing devices, which cannot be said about the use of differential null indicators, 


Here we consider the stability of one type of group-one circuit, namely a four-arm bridge with two identical 
servodrives which actuate the resistive and reactive components in one arm (Fig, 1). 


Krasovskii's complex coordinates method [7] can be used in the analysis. 

We deal here with the effect of the static nonlinearity and of the balancing circuit coupling on the 
stability, 
2, The Stability of a Linearized System 


In drawing up the equations of motion for the circuit of Fig, 1 we neglect the transients in the bridge 
proper (see Appendix) and use the measurement diagonal voltage in complex form: 


At small deviations from balance the static characteristics are arcs of circles and can be replaced by 
their tangents at the balance point, Let subscript 0 denote values at balance; then 


AU = AZ, + (=) AZ, = —k, (AZ, — AZ,), (2) 
0 0 
where 
k = k,e-?% Jiohts e—i% 
V (Rp + 
U X10 (3) 
AZ, = AZ, Fe is the referred value, 


The primes are in future dropped, but we bear in mind that AZ, (or ARy and AX,) appearing later are 
actually referred quantities, 


If the reference voltage of phase-sensitive detector I differs in phase by an angle gy from the current in 
branch one (iyo), and that of II by an angle @ = y+——, the output voltages €, and ¢€3 from the detectors when 


these operate as switches are,respectively, 


e, = [AU — kykgRe AZ,)], 


(4) 
’ = ky Im [AV = — Im — AZ,)], 
where kg is the gain of P and of the detectors, 
Since AZ; = AR + JAX, and AZ, = AR, + jAX, (4) becomes 
[cos (9 _ (AR, — AR;) sin (9 (AX,— AX;)}, (4a) 
= — k,k, [sin (p — a) (Alt, — AR,) + cos — a) (AX, — AX,)}. 
Adding to (4) or (4a) the equation for the servomotors plus amplifiers and balancing elements 
(Tp +1) pAR, = kskye,, (Tp +1) pAX, = (5) 


where p — , kg is the amplifier gain, k, is the gain of servomotor plus reducing gear and balancing element 


and T the servomotor time-constant, we get a system of equations describing the action in the linear approxi- 
mation, The voltage amplifiers and detectors can in most cases be taken as lagless, 


Figure 2 shows the equivalent structural diagram of an automatic bridge compiled from (4a) and (5), 
which shows that the bridge is a two-channel servo with direct antisymmetric dc feedbacks with relation to the 
input parameters (AR, and AX;) and to the controlled parameters (AR, and AX,). A special feature is that the 
coupling constants k, = sin (g— a») can be altered via the phases of the reference voltages, 


“4 taal 


as, ah, 
&Y 
Ky 
\ 
Fig, 2 


We now examine how the coupling constants affect the properties of the system, To produce the necessary 
sensitivity in measuring the input parameters we must have a gain given by 


U 
K = kykgks cos (p — a.) = (6) 


where U, is the servomotor starting voltage and § the minimum change in the input to which the system must 
react, 


Equation (6) shows that as the coupling constants (i,e., |g — a |) increase, kykgks must be increased to 
keep K constant, 


Hence the open-loop gain Kp = kyksksky cos ( ¢ — a) must be independent of |y— ay|. 
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The stability condition implies that we must have . 
Ky<Kper, (7) 


where Kp. cr 28 the critical open-loop gain, 
We use the complex coordinates method [7] to determine how K,, ., depends on the coupling constants, 
Combining (4) and (5), and applying the Laplace transform (initial conditions zero) we have 
(TS + 1) SAR, (S) = —k, Re {{AZ,(S) — AZ, 
(TS + 1) SAX, (8) = — k, Im {[AZ, (S) — AZ, 8) 
where kg = kykgkglky, and AR, (S), AX4(S), A24(S), and A2,(S) are the Laplace transforms of ARy, 424, 
and AZ. Multiplying the second equation in (8) by j and adding to the first, and using the fact that AR, (S)+ 


+ jAX,(S) = AZ, (S), we get the equation for the system in complex form: 
[(7S + 1) S + AZ, (S) = AZ, (S). (9) 


The closed-loop transfer function is 
Ww, (S) 
; (10) 


where 
Wy (S) (11) 


is the open-loop transfer function, 
For stability 1 + W,, (S) must have no zero in the right half-plane of S, and so Wp (jw) must not encompass 


the point [— 1, j0] when w varies between — @ and +, 


Equation (11) implies that the circuit coupling amounts to rotating the function We (jw) = 
(the complex transfer function for the separate balancing circuits, with = a) through an angle ( g — a») de- 


(Tjw +1) jw 


termined by the degree of coupling, 


H 
0 
‘ 
au 


Fig. 3 


For this W, (jw), and also for other types found in automatic bridges, rotation about the origin means 


that the hodograph of W, (jw) approaches the critical point{— 1, jo), i.e,, that Kyo; is reduced (Fig, 3), 
Hence Kp cr decreases as ke increases, while to keep the sensitivity and response time (which depend 


on ky) constant, Kp must be kept constant, 
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This leads directly to the conclusion that the balancing circuits should be uncoupled, 


As the circuit of Fig. 1 has nonlinear static characteristics the balancing circuits can only be uncoupled 
(at small deviations) by appropriate reference voltage phase selection at one point in the working range, At 
other points a will differ from ag = ¢. 


The point where gy = a is selected within the working range such that the coupling at all other points 
is the least possible, 


We must point out here that these bridges can be designed to have balancing circuits uncoupled over the 
entire working range, This can occur if we connect the balancing elements as in Fig, 4[8], i.e,, in such a way 
that the current in the arm remains fixed, Here the static characteristics (balancing lines) are straight, as the 

topographical diagram of the adb branch (Fig. 5) shows; un- 
coupling at one point implies uncoupling at all, 


An additional advantage of the circuit of Fig. 4 over 
that of Fig. 1 is that capacitative and inductive reactances can 
both be measured without switching, since the potential point 
corresponding to d can move into either the upper or lower 
half-planes in Fig, 5 as the balancing elements are adjusted, 


A fault often found in the circuit of Fig, 4 is that the 
readings are not separate, and are frequency-dependent, How- 
ever, these features are not essential in certain applications, 

An example is the checking of dimensions and conductivity of 
parts using eddy currents with an inductive transducer (9), Here 
the checking method is frequency-dependent, The frequency 
must be known, and the two measured parameters, e.g, diameter 
and conductivity, are both related to the resistive and reactive 
components of the transducer impedance, so the two parameters 
cannot be read off separately on a normal bridge either. 


Fig. 5. Topographical diagram for the - 
adb branch, 1) Balance line on adjusting 
the resistive arm; §R=var, 6X =const; 
2) balance line on adjusting the reactive 
arm, §6X=var, §6R=const, The points 
4, e, d, g, u, and b in the diagram corre- Characteristic Nonlinearity 
spond to the points of a, e, d, g, u, and b 
in the bridge. 


3, Stability Allowing For Static 


Consider the behavior when AZ, is large but does not fall 
outside the range in which coordinate bounding can be neglected. 


Putting in (1) AZ, = AZ and using the balanced bridge equation 


= 


we get 
—Z,AZ 
Z3 + Zo +42, (12) 
From (4), (5), and (12) we have 


(Tp +1) pAR, =k, Re 


— = F,(AR,, AX,), 
Rs + + AZ, | 
AZ, 


Rs + Zy + AZ, 


(13) 


(Tp +1) pAX, = ke Im| F, (AR,, AX,), 


where , 
F3(AR,, AX,) = hi (AR,, AX,) cos? — he (AR,, AX,4) sin 


Fy(AR,, AX4) =f, (ARy, AX,) sine + fy (ARy, AX,) 9, 
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(12s + + AR, 4- (Xan + AX, 
(Rs + Ray + AXg— (Xu + AX,) AR, 


ks kikok sks, ky 


Putting (13) as a normal system of equations, with notation as below, we have 


AR, = AX, = Zo, 


dz 

= (15) 


Our interest is directed to the behavior with the initial conditions 
= = AX$, 28 = 2° =0. (16) 


We restrict consideration to cases where T can be assumed sufficiently small, as is true of some designs 
of automatic bridges, 


Then, using the results of [10], we can replace (15) by the degenerate system 


d. 

d (17) 
=0 


with the initial conditions 
= 22 = AX}. 


In fact the right halves of (15) are bounded and continuous in the xj(i=1,..., 4), and define the 

functions Xs = Fs (x3, X4 = F4(xq, also continuous in the xj, The equations = — xj +F;=0 (i=1, 2) 

define the stable root (since = = — 1<0, i=1,2) and these points in the region covered by (16) are points 

i 

internal to the region of influence of the stable root F’; = 0 (i=1,2), Hence, from [10] (p. 202), when T—> 0 
the solution to (15) subject to (16) tends atthe limit to the solution to the degenerate system (17) subject to (18), 
this limit being stable relative to the initial conditions of (16). 


(17) can be solved, using Erugin’s method [11], by introducing the complex function AZ of the real 
variable t: 


AZ, (1) = AR, (t) + FAX, (0). (19) 


Multiplying the second equation of (13) by j and combining with the first, allowing also for (19), and 
then putting T =0, we get the trace of (17) as a single complex equation: 


dAZ, _ K AZ, 
at Rs + Zig + AZ, 


(20) 


Since the RHS of (20) is an analytic function of AZ, for all possible values of the latter, integrating (20) 
and using (18), we get a solution to (20): 


Rs + Zo AZ, 1 3 40 
Az? + ke’? (21) 


where AZ? = + jAX%. 
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Putting (21) in exponential form we have 


ws AZ, Ze (29) 
where 
— Rat Lu 
(22a) 


(22) shows that the sole balance point when Re + =a<0 is AZ, =0 (AR, =0, AX,=0). To show this 


it is sufficient to note that at t=oo and a <0 (22) is satisfied only when |AZ, | =0. For all AZ, | # ow this 
follows directly from the relationships for the moduli of (22) 


cos (y— —R 


| 


where 


k 
Y = Re + cos (p — a»). 


To prove that AZ = 0 is not also a root of (22) when t = w (a < 0), consider the limit 


_ AZ, 
|AZ,|-+00 |AZ,|->00 
The modulus of this limit is 
cos 
lim eRst Ze lim Rat A (25) 
|AZ,|->00 


When |AZ |—> 0, cos (y —a) > 0, This is readily seen by considering Z, in the complex plane 


(Fig. 6), Thus when |AZ,| —> o the RHS of (22) tends (in modulus) to infinity, i.e., | AZ, | =o is nota 
root of (22), 
Hence, when a < 0, any initial deviation A y eventually 

iM Jar : (4) results in a return to balance, The condition a < 0, which amounts 

y- to the condition 

cos (p — >0 (26) 

; H [since k < 0; see (14)], is therefore the necessary and sufficient con- 
— aR a dition for stability “in the large” for this type of bridge, 
fu Equations (13) were compiled for the balance point at which 
| 
R — 
‘ arc tan 
Fig. 6 


If the bridge is to be stable for all possible Ry and X, we must 
have cos( g— a) > 0, where the range in a is determined by the ranges in Ry and X,, Since in any bridge 


0< |a|< ‘ we can always ensure stability for all possible Ry and X, by suitable choice of ¢, 
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Condition (26) for the stability “in the large” is the same as the condition for the linearized system, 
Hence the static nonlinearity does not here lessen the stability. 


Equation (21) enables us to determine either the maximum time to reach balance for a given g, or the 
g for which the maximum time to reach balance is least, 


APPENDIX 


Neglect of Bridge Circuit Processes 


Each balance circuit in the bridge of Fig. 1 contains modulated sections (bridge M, preamplifier P, and 
detector N), 


To be able to use the methods employed with ordinary circuits (unmodulated) we must first establish the 
equivalent transfer function in terms of the envelope for the modulated section, This section has two inputs 
and a lag in the modulator, 


¥ 
Fig. 7 


Consider one bridge branch (Fig. 7), which is a modulator with a lag, We suppose the balancing element 
to be an inductance, The deductions below are equally correct if the reactive balance element is a capacitance, 


The instantaneous branch parameters are given by 


d 
dt (Lyx) + + Rei = Us (27) 


Linearizing (27) about the balance point (Ry = Ryo, Ly = L4o, ig = igo) and supposing that 


ig ¢, (28) 
we get the branch response equation as 
(tp +1) AU, =—k(AR, += pAX, sine t—kAX, cose t, (29) 
where ‘ 
t= = k = Ih Re’ AU, = Ai, Re. 


From (29) the branch can be considered as a modulator, whose output is the modulated carrier 
(t) + sinw t+ AX, (t)cosw t= f, (t)sinw t+ f, (t) cos @ 
in series with a linear link of transfer function 
(30) 
Hence the problem amounts to that of passing a modulated carrier Uj, = fy(t) sin wt + f(t) cos wt through 
a linear circuit, The complex transfer function method given in [12] is convenient to use here, From [12], we 


introduce two independent operators j and i such that i? = j*=— 1, but ij=jiz#— 1, Operator i refers to 
the signal envelope and j to the carrier. 
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The output signal from the linear link, in complex form, is 


Uy = hi (t) + /fa(t), 


or after a Laplace transform 


(S) = Fi (S) + jF2(S), 


where 


1 
Fy (S) = (t)] = (S) + S4X1(S),  F2(S) = L [fa (t)] = AX, (5). 


The variable transform S (envelope) is only complex in the i plane, 


The input signal to the linear link of transfer function K(S) is defined by 


out (S) = AU, (S) = K (S + jo )U,(S) = 
1 
= K(S+jo ) (S) + SAX, 


(33) plus (30) define both complex transfer functions for the envelope for the branch 


AU; (S) —k —k[(tS + 1) — | 


Wal5) = RRS) FIO +1 +o)?’ 


AU, (S) 


1 
+ i) (es + 1)+to + i 


Wx (S)= TS Fie) +1 (tS + 1)? + (to 


The equation for this branch will be 


AU, (S) = Wp (S) AR, (S) + Wy (S) AX, (5). 


For the other branch we have an analogous equation: 


AU (S) = Wp (S) ARg(S) + Wy (S) AX4(S), 


where AR, (S) and AX,(S) are referred quantities (see Section 2), 


where 
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The voltage on the measurement detector diagonal will be 


AU (S) = AU2 (S) — AV,(S) = Wp (S) [A Rg (S) — (S)[ + Wy (S) [AX (S)— AX, (5). 


The output voltages from the detectors are given by ~ 


(S) = Rej [AU (S) = 
= Wp (5S) [ARs (S) — AR(S)] + Wy (S) ([AX4(S)— 4X, (S)], 


(S) 


= Re, [AU (S) kee?®] = 


= Wp (S){ARa(S) — (S)] + Wx (S) [AXa(S)— 4X, (S)], 


WAR) = Res [Wp (S) = 


(S) = Re j [Wy (S) = 


— kkg[(tS + 1)cose +t sin 9] 
(tS + 1)* + (tw 


— kk S (tS + 1)+ to ] 08 —sin 
(tS + 1)? + (tw )* 


(31) 


(32) 


(33) 


(34) 


(34a) 


(35) 


(36) 


(37) 


(38) 


The expressions for W", (S) and W"y(S) are obtained by replacing g by @ in (38), 


The functions W'S), W'y(S), W",(S), and W"\(S) are equivalent transfer functions for the measuring part of the 
bridge, which is the modulated section, 


Supplementing (37) with the servodrive equations 
AR,(S) = Ki (S) (5), AX4(S) = K2(S) €2(S), 


we get a system of equations describing the automatic bridge: 


Fens + Ap (s)| AR, (S) + (S) AX4(S) (S) (S) (5) 4X1 (5), 
39 
(S) ARa(S) + | + Aig (S)] = (S) (S) + AY (S) (5), 
where 
D(S) = + 275 + 1+ Q?, 
Ap (S) = — kks [(tS + 1) cos + Qsin 9], (40) 


Ay (S) = —kk, (tS + 1) + cos e—sing| (Q = Tw ). 


The expressions for A"p(S) and A"y(S) are derived from those for A'p(S) and A'y (S) by replacing @ 
by @. 


The characteristic equation for the system is put as 
D*(S) D(S) ° D(S) ’ 
Ki(5) Ki) = Ar (5) + Ax(S) Ki(s) Ap(S) Ay (S)— Ap (S) Ax (S) = 9. (41) 


In most cases of designing ac bridges Ryo + Rg and X,q are chosen of the same order, because comparable 
sensitivities in both quantities measured are required [3], Hence 


is of order & « H, i.e.,even at 50 cps it is of order 0.01, 
wW 


We now examine how the stability can be judged from the degenerate characteristic equation obtained 
from (41) by putting r =0, 


Here it is necessary [13] to show that an accessory equation is stable, 


The structure of A (S) and D(S), which are polynomials of identical degrees in S and r, show that the 
characteristic equation is of the form 


[MogS™ + My + ... Moy) + + + 


(42) 
+ M,,SN—"-1 + = 0. 
Hence, the appropriate accessory equation must be a first-order one, i.e., 
M9” + + +M,, = 0. (43) 
Suppose 
P,(S) 
KS) = SQ, (sy = » 


where P,(S) and P{S) are polynomials of degrees J, andJ,, and Q,(S) and Q,(S) are polynomials of degrees 
my and mg, with <my and mg, 
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Oy Rig + Re Oy Rig + Rg 


We then see that the coefficients to the highest powers of S for each power of r in the characteristic 
equation, i,€., Moo» Myo, . «+ » Mno in (43), coincide with those for the corresponding degrees of rS in DXS), 
Hence, the accessory equation is obtained by equating D*(S) (with r3 replaced by q) to zero: 


D? (q) = (q? + 29 +14 = 0. 


This equation always has a stable solution, 


Hence the stability “in the small" for automatic bridges with astatic balancing circuits in which Xy and 
Ryp + Rg are of the same order can be judged from the degenerate characteristic equation obtained from (41) by 
putting r =0, i.e., we can neglect transients in the bridge circuit, 
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FREQUENCY METHODS OF REMOTE CONTROL FOR DISTRIBUTED OBJECTS 


V. A. Il'in and K, P, Kurdiukov 


(Moscow) 


The analysis leads to the selection and development of a frequency method of selecting and controlling 
distributed objects using line control channels, The design principles for remote control devices using frequency 
relays with series-resonant circuits are considered, 


The methods and techniques of remote control so far developed for power systems, aircraft, railroads, etc, 
have mainly been for lumped objects. 


As better process control methods are introduced the need arises to deal with many small control points 
or isolated objects distributed over areas or along lines (e.g. in oil- and gas- production, in irrigation systems, 
large factories, mines, pipelines, etc.), 


The number of such points dealt with from one center may be tens or hundreds, while the total number 
of control or telemetering signals is small, seldom more than 3-5, 


As the lines commonly used in such systems are relatively lightly loaded, and the signal to noise ratio is 
high, efficiency of information transfer and electrical noise effects are not commonly of major importance, The 
main problems in such systems are reliability in apparatus and lines, plus simplicity, convenience, etc, 


The apparatus at the central station is usually under continuous observation, while the remote objects 
(outstations) are often in inaccessible places or in closed boxes in the open, and are left unattended for long 
periods, Hence, these latter objects must be much the more reliable, and the reliability of the system as a 
whole is largely determined by that of the outstations, 


Various selection methods are considered here; the frequency method of selection and control is chosen, 
The design principles for remote control devices using frequency relays with series- resonant circuits are con- 
sidered, The choice of line structure for distributed objects has already been dealt with (1). 


Choice of Selection Method 


The reliability at the controlled objécts mainly depends on the properties of the electromagnetic relays, 
tubes, gas discharge devices, contactors, etc., as these are relatively liable to break down, 


In the devices to be considered no tubes or gas devices are used, so the most unreliable elements can be 
considered to be the relays. Hence, in choosing methods the number of relays involved should be of primary 
importance (other things being equal) neglecting for the moment the numbers and loadings of the contacts, 


The uses of contactless relays for this purpose have as yet been little studied, so these are not considered 
here, Their use undoubtedly ought to be examined and extended, 


In choosing a selection method one must also ensure identity in the control equipment at the outstations, 
since this facilitates manufacture and servicing, 


Remote control and telemetering devices for lumped objects (power grids, railroads, underground railways, 
etc.) often use time-division with the signals (e.g. in devices of VRT, RVS, KRV, DVK, etc. type which use time 
selection markers [1]), Attempts have been made to use the same method with distributed objects, Here each 


167 


effector point must control a rack of equipment with a distributor adequate to carry the full load of the system, 
although each such point may contain only one or two of the tens or hundreds of such objects in the system. 
Complex selection devices are then needed, 


Let us determine the number of distribution and selection relays per object using time division methods 
with varying selection systems, We only consider a single line with a total of N objects connected to it (2). 


Using time division methods the number of distribution and selection relays per object, k, depends on 
the selection method; in general k= m+n, where n is the total number of pulses in the code and m the number 
of selecting pulses, 


If we neglect economies made by using repeated relay operations during a single transmission, rectifiers, 
etc,, the number of distribution and selection relays per object can be described by the following formulae, 


1, Distributive nongrouped method (N =n, m = 1) 
k=N+1, 


2, Distributive grouped method with q degrees of selection and total number of objects per group N vq 
and m=q 


q -—- 
k =q(VN +}). 
3. Combination selection, with k=n+m for one set; then n, m and N are related by 


N n! 


(a — m)\m! f 


4, Combination selection, with for all sets 


n 
mc™ 

log 2 


5. Distributive grouped method with distributors located at the objects, and with a number of pulses in 
the code equal to the number of objects; the line sections are connected to the group point one after the other 
when the code is being transmitted [3], and each object has a different number of code pulses: 


K 
26 7 7 
2% 
22 / 
20 
8 / 
6 
6 
4 Pad 
5 
2 
67890 20 «30 50 00 200N 


Fig. 1, Relation of the number of distribution and selection 
relays per object to the total number of objects with the various 
selection methods, 


. 


n 
k= 
Figure 1 shows k= g(N) for the various selection methods, 


Figure 1 also shows the number of distribution and selection relays per object for the frequency method 
(curve 6, k = 1), each object having a distinct frequency, The relays are here the effectors, and play no part 
in selection, 


Figure 1 shows that the number of distribution and selection relays per object can be minimal for the 
frequency method with distributed distributors [3]. 


All other methods require many more relays per object and (more important) the number of such relays 
depends on N, 


Hence, using distributed distributors and frequency methods we can produce standard apparatus suitable 
for varying numbers of objects, Both methods have the substantial disadvantage that it is difficult to check the 
situation at the objects, particularly if the lines are long and the objects many, 


The first method has the following disadvantages, 


1, The lines and system generally are less reliable because they contain many contacts (at least N) 
connected in series, A fault in one relay puts the whole system out of commission. The structural reliability 
is low, 


Line point Line point 
N2 
GP 


Fig. 2, Line point circuits using frequency selection and polarized 
control pulses, 


2, It is difficult to prevent noise distorting the code, if a noise pulse transforms a counting pulse to a 
selection pulse, If all relays operate equally rapidly the selection time is some N times larger than in the 
frequency method, The probability of interference from noise is correspondingly increased, 


3, It is difficult to ensure independent calls and telephone conversations using a single line, although 
this facility is in many cases necessary in such telemechanics systems, 


These disadvantages restrict the use of the system, The frequency method also has disadvantages, the 
main one being that interference can produce distortions after selection (during the control or telemetering 
operations). This effect can be reduced to negligible proportions by using coding methods, For example, if 
the object is selected by frequency and polarized pulses used for control (e.g. for switching on the telemeter 
transducer) the relay can be arranged to be cut in by the pulses of one polarity only, and to be held on by 
current of reverse polarity, Figure 2 shows this system for a one- wire circuit, Here GP is the group point, the 
Fy (K =1, 2) frejuency selection relays, and the Yy are control relays slugged on release, If a noise frequency 
occurs during the control time (after selection) and causes the frequency relay in some other object to operate 
(false selection), no undesirable or dangerous consequences arise, since a direct current is supplied to the line 
during the control period of sign opposed to that on which the control relay will operate, 


A notable feature of frequency methods is that they are very flexible and that the various objects are 
virtually independent, If adequate margin is available in the frequency band, new objects can be connected 
without disturbing the old, 
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The frequency method is therefore often used with distributed objects, particularly if the line lengths 
are not more than a few tens of kilometers, or the number of elements per Wire not over some dozens, The 
branched line structure being used make the use of this method more suitable, 


Frequency Relays With Series Resonant Circuits 


The frequency relay, consisting of relay plus frequency selection device, is the main component in 
frequency remote control and telemetering circuits, 


Line point Line point 


GP - ----- 


=> 
= 


Fig. 3. Line points with series-resonant circuits and neutral 
relays. 


Mechanical resonant systems (tuning forks, reeds, strings) electromagnetically driven have been much 
used for frequency selection; these either cause leaf spring contacts to close directly, or else operate polarized 
relays, magnetic amplifiers or transistors via rectifiers connected to their outputs, These systems have narrow 
response bands and are relatively stable [4-6], A disadvantage is that the efficiencies are low and frequency re- 
setting is difficult. Individual adjustment is required, the inventory has to be large, and the cost is high, as mass- 
production techniques cannot be used, Much better results are had by using series-resonant circuits with Q= 


r 
K 
In the audio and subaudio ranges inductors built on standard toroid cores are best used, The toroid acts 
simultaneously as a stepdown autotransformer, with taps for selecting the operating voltage of the relay connected 
to it via a bridge rectifier (Fig. 3), 


Figure 4 shows the experimental Q = ¢ (f) curves (f = frequency) for coils wound with various wires on 
AlSiFe TCh-60 cores, The use of PEV wire instead of PELShO increases the Q by 25 to 40%, The impedance 
at resonance is usually made several times that of the line, The relation of load and resonant impedances 
is determined by the permissible reduction in the Q of the unloaded circuit, If the frequency is low and 


< 

where the subscript 2's relate to the line, the simplified equivalent circuit for one such system connected to 
the line will be as in Fig. 5, Here rj is the internal resistance of the source, Ix =Ict ty is the resistance of 
the circuit at resonance andr the load resistance as seen by the circuit. 


For an above-ground two- wire line 10 km long consisting of steel wire 4 mm in diameter operating at 
200 cps we have 
ob, = 1256 x 9x 10% = 110 ohm, r, = 140 ohm, 


@C, {256x102 = 16 kohm 
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The circuit can be designed to this line from 


oL,, 
Q 
150 
2 


Y | 
25 7 
| 
20. 40 50 100 200 400 1000 f, cps 
Fig. 4. Variation of Q with f for the toroids, Data on 
Toroids : 
Curve No.| No, of Type of| Wire dia- | Inductance, | Resistance, 
turns wire meter, mm | henries, ohms 
1 6500 PEV 0.31 6,15 125 
2 6500 PEV 0.33 4.12 100 
3 5500 PELShO 0.31 3.86 120 
4 4300 PELShO 0,31 2.47 75 


The efficiency at resonance with wl « ty for very low frequencies is given by 


r 


nth +? 


{= 


In some cases it is desirable to use at very low frequencies the coefficients 


r r; r 
min 

ky = r ’ ky = r ’ ks r 
K min K min « min 


where Tt) min is the line resistance between the group point 


‘ f and the nearest object, and to assign the frequencies in such 
a way that the frequency increases with distance, In this 
r way the reduction in efficiency at the higher frequencies 
can to some extent be compensated, 
i i the effic nd can 
Fig. 5. Simplified equivalent circuit for Being given ky, ke, and ke, the efficiency and Ly 
be determined from theQ= g(f), 
one system connected to the line. 
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The frequency relay resonance parameters determine the over-all Q: 


It, at the working frequencies is found from Q = g(f), Then Q, = g (f) curves are drawn up to determine 
the minimum spacing between frequencies, allowing for line parameters and supply voltage. 
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Fig. 6, Line point with a resonant circuit for 


Fig. 7. Line point with a resonant circuit for 
selecting combinations of two frequencies, 


selecting combinations of three frequencies, 


Higher- resistance circuits are used for the longer lines, or with more numerous objects. 


A Frequency-Combination Selection Method With One Relay Per Object 


Electrically resonant relays enable one to use frequency-combination selection with one relay per object 
only; this method extends the use of the frequency method, When this new method is used with a group point 
in the line, two frequencies fy and f, are transmitted together or successively. f, is supplied briefly, and f, until 
the end of the operation, When f, arrives the relay, tuned to f, in the quiescent state, operates and switches the 
circuit over to fz, by replacing Cy by Cg, as in Fig. 6. f, holds the relay on throughout the operation, 


GP 


Line 


R Selection 


relay 
Fig. 8. Simplified telemechanics circuit using 


frequency selection and polarized signals. 


Fig. 9. Simplified telemechanics circuit using 


frequency selection and frequency interrogation 
The number of objects N (code formula) in this 
‘ of the position from the control, 
system will be nyng, where ny is the number of f,'s and , 
: C,, is a preset capacitor to tune to the working 
nz the number of f,'s, P 


frequency, C, an additional capacitor, 
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To increase the reliability the relay must be a slugged one, or the integrating capacitor C must be used, 


The discharge time-constant for C must be such that during the transfer period the voltage on C does not drop 
below the hold-on value, 


This method can be used with three or more frequencies also, If three are used all the apparatus needed 
at the object is one inductor plus two electromagnetic relays 1 and 2, as shown in Fig. 7. 


Frequency and polarized signals can be used simultaneously for selection. In this way circuits normally 
disconnected from the line can be used, as in Fig. 8. A brief dc pulse and the relevant frequency are dispatched 
together or successively from the group point. 


The dc pulse passes through one winding of the relays R in all selectors, and causes them to operate, A 


relay can only hold on if its resonant frequency is also supplied, This greatly reduces the circuit loading of the 
line, 


Remote Checking of Object Positions 


Frequency relays using electrical resonance provide simpler means of checking object positions, 


Relay volts 
ff 
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Fig. 10, Resonance curves for different toroidal coil windings connected to 
the load (relay), The numbers are those of the turns, 


173 


With short lines (up to 5-10 km) and ry > rp the object position is best checked from the group point 
using the change in amplitude of the frequency used for selection, Figure 9 shows a circuit for this purpose, 


Relay volts 
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Fig. 11, Resonance curves for different toroidal coil windings 
connected to the load (relay). 
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TABLE 
No. of Circuit No, of turns | resonant Line resistance to point, ohms 
curve on 
group seconda frequency, 
uf ry full broken dot and 
- curve curve dash curve 
1 500 25 0 
2 4.5 500 34 0 — _ 
3 2.8 500 44 0 _— _— 
4 2.27 370 54 0 11 15 
5 1.32 370 71 0 14 24 
6 0.74 275 92 0 17 29 
9 0.24 200 157 0 27 50 
10 0.17 200 182 0 30 57 
11 0.138 200 210 0 33 64 
12 0.106 140 240 0 36 71 
13 0.08 140 269 0 39 78 
14 0.063 140 300 0 42 85 
15 0.05 140 335 0 45 92 
16 0.039 140 373 0 48 _— 
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Relay volts 


Fig. 12, Resonance curves for series circuits evenly spaced along the line, 


0 10 2% 30 40 50 60 70 30 100 120 130 150 160 1 180 190 200 20 220 280 260 250 260 210 280 290 300 HO 320 10 380 310 cps 


The contact on the object switches an additional 


capacitor into or out of the circuit at the extreme positions, 
and thereby alters the resonant frequency, One frequency 

is the selection frequency, while the other can, for instance, 
be common to all objects, or all of a group, If the object 
is, for instance, switched on, the extra capacitor is cut out 
and the circuit has its proper frequency fy, The resistance 
at resonance is reduced, and so the check relay PK at the 
central panel operates when fy is transmitted, If the object 
is in the other position the extra capacitor is cut in, and 

the other frequency (e.g. f,) is adopted; and so the check 
relay at the central panel does not operate when fx is 
transmitted, since the impedance at this frequency is higher, 


By transmitting a brief pulse of frequency f, and 


also fy, the object can be cut in, and the switch-on checked, 
With lines over 5-10 km long this distant check is best done 
using the contact on the object to operate a local generator 
operating at the frequency proper to the point, This generator 
can be built round a transistor fed from a local battery or 
other source, The check can also be performed by using 
damped oscillations at the resonant frequency excited by a 
pulse sent from the control point. 


EXPERIMENTAL RESULTS 


Figure 10 shows resonance curves for a circuit as of 


Fig. 3. The relay voltage was measured at the output of 

a bridge rectifier connected to a RKN relay of resistance 
350 ohms, The line voltage was constant at 20 volts, The 
toroid was wound with 0,31 mm diameter PELShO wire on 
an AlSiFe TCh-60 core, In all 5500 turns were wound, 
tapped at 750, 500, 370, 275, 200, and 140 turns, C = 

= 0.75 yf, Curves were then taken for the load placed 
across various taps, Figure 11 shows the same curves on 
the relative scale (Relay volts/Max, relay volts) = ¢ (f), 


Figures 10 and 11 show that as the resistance in- 


troduced by the load increases (with the number of turns) 

the resonance becomes flatter, The ferroresonance effect 

is also reduced, so the curves become less asymmetric. 

By using a working frequency somewhat below the resonant 
one, e.g. 90 cps, the ferroresonance effect can be neglected; 
this effect reduces the certainty of relay operation if the 
voltage alters, 


Figures 10 and 11 show that as the damping increases 
(including effects due to increased line resistance) the 
method becomes of more limited value, 


Figure 12 shows resonance curves for series circuits 
evenly spaced along the line, 


The line and circuit parameters for each of the curves 
of Fig. 12 are given in the Table. In all cases the induct- 
ances were equal; the data have been given above, 
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The maximum resistance of the line from the generator to a circuit was 92 ohms, which corresponds to 
about 8 km of 4 mm diameter steel wire, 


The demands as to the generator at the central point are comparatively easily satisfied with these circuits. 
Variations of 1-2% in frequency and 5-10% in output voltage are permissible, Harmonics are permissible up to 
10-12%, 


SUMMARY 


1, The frequency method of selection results in least relays per object, which greatly reduces the total 
number of relays in systems with distributed objects. The remote control apparatus at the object can then be 
very simple, reliable and standard, in spite of variations in system size, This makes the method very suitable 
for systems with distributed objects, 


2. Interference can cause a second object to be selected and controlled after a first has been selected 
in the frequency method, This basic disadvantage can be largely eliminated by using combined selection 
systems, e.g. frequency plus polarity or frequency combinations, 


3, The LC serfes-resonant circuits described in the paper can often be used to advantage in systems with 
distributed objects as the frequency selection devices, 


4, The frequency relay described can be used in frequency combination selection methods, One relay 
per object is then required, 


5. The amplitude-frequency of checking the object position described above enables one to simplify 
the signalling system for lines up to 5-10 km long. 
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LETTERS TO THE EDITOR 


A NONLINEAR SEMICONDUCTOR OF RESISTANCE SENSITIVE TO 
MAGNETIC FIELDS 


G. E, Pilus and O, V. Sorokin 


(Leningrad) 


Several theoretical and experimental papers have been published on resistance changes in semiconductors 
in magnetic fields [1-6], It has been shown that in semiconductors with carriers of two signs change in carrier 
concentration occurs on passing a current; this is additional to the normal increase in resistance (pioportional 
to the square of the field) caused by the change in carrier mobility, If the current in the specimen and the ex- 
ternal field are directed as shown in the figure and the specimen ts a thin wafer of thickness d, then the total 
change in resistance will under certain conditions [1-6] be given by 


ws, (1) 


where fp is the resistance at zero electric field, b = H p/H n Y = ap/ap (where Up and Hp: are the electron and 
hole drift mobilities and ap and ap are the ratios of Hall to drift mobilities for electrons and holes), U is the 
voltage applied, H the magnetic field strength, d the thickness, 7 the length, T the absolute temperature, k 
Boltzmann's constant, ande_ the absolute electronic charge, 


Consider, for example, the characteristics of a nonlinear germanium resistor of inherent conductivity 


n=p=2.3x10%cnr*, 4, =3700cm*/sec, Hp=1700 cm/sec, a,—1, y= 4.2, 
with a carrier diffusion length ,9 =0.2 cm and specific resistance p =50 ohm-cm, We putd=0,15 cm, | = 
=1,.5 cm, h =0,.5 cm, 


y For room temperature (1) can be put as 


Ar = 0.083 x 10-8, UH, (2) 


Ye Ws where Ar and fp are in ohms, U is in volts and H in oersted, 


If, for instance, a permanent magnet with H = 1200 is used, 


then (2) gives 


r=ro(1+0.1U). (3) 


A useful feature of this device is that the sign of the resistance change depends on the field directions, 
Since the resistors are small, a block of them, appropriately orientated and connected, can be placed together 
in the working gap of a single magnet. 


The nonlinearity can be increased by replacing the permanent magnet by a magnet fed from the same 
source as the resistors, In this case the variation of Ar with H is also utilized, The resistance can also be ad- 
justed by varying the magnetic field independently. 
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CHRONICLE 


THE INTERNATIONAL FEDERATION OF AUTOMATIC CONTROL (IFAC) 


The Two First Periods of Development of IFAC 


The International Congress of automatic control was held in Heidelberg (W, Germany) in September, 1956, 


The congress was organized by the Control Group of the Union of German Engineers (VDI/VDE, Dusseldorf); 
the president was Dr, O, Grebe, Some 1000 representatives of 18 major states attended, including many famous 
engineers and scientists, 


The soviet delegation, composed of A, M, Letov (Doctor of Physico- Mathematical Sciences), la, Z, Tsypkin 
(Doctor of Technical Sciences), V, N, Naumov (Candidate in Technical Sciences) and V, A, Racheev (Engineer) 
presented papers, 


On September 27, Dr, Grebe called together representatives of the countries participating to discuss the 
proposal from Professor Oldenburger (U,S,A,) that an International Federation of Automatic Control be organized, 
with aims: 


1) to facilitate the exchange of scientific information on automatic control between member countries, 
and 


2) to organize international congresses on automatic control every four years, 


The Soviet delegation gave its full support to this proposal, as extending contacts between scientists and 
engineers in different countries, 


The aims of the federation were generally agreed, and the congress elected a Preparatory Committee 
composed of Professor V. Broide (France), Dr, Wellburn (Great Britain, president), Dr, O, Grebe (W, Germany), 
P, Nowatski (Poland), Professor Oldenburger (U.S,A.) and Dr, A, M, Letov (USSR), 


The main task of the preparatory committee was to organize IFAC, 
The work was divided into two sections: 
1) establishing the views of countries on the federation; 


2) developing proposals for the structure of the federation, for its statutes, and for calling the first General 
Assembly of IFAC, 


In the USSR the first part was undertaken by an interim national committee set up by the Academy of 
Sciences, 


The Preparatory Committee drew up a proposed charter for the federation, defining its aims, methods and 
means, membership, structure, rules and obligations on members, and also defining the functions of the General 
Assembly, 


These statutes were drawn up using two alternative sets proposed respectively by Professors Broide (France) 
and Oldenburger (U,S,A,), and were adopted at a special session of the preparatory committee on April 25-27, 
1957, in Dusseldorf, The resolution calling for the first General Assembly, and drafts of many resolutions to be 
presented at it, were also adopted, 


The second meeting was set for September 9, 1957, in Paris, with the first General Assembly to follow on 
September 10-12 at the same place, 


The section on membership occasioned the most discussion, The Soviet delegation considered that any 
country which expressed a desire to join IFAC should be admissible, Another view was that the membership 
question should be decided by vote at the next General Assembly, This was entered in the statutes, 


In addition the committee agreed that states which notified their acceptance of the statutes of IFAC 
in writing prior to September 9, 1957 (i.e, prior to the first General Assembly) should be members, 


Meeting of the Preparatory Committee September 9 and 10, 1957, 


This meeting was held with Professor Broide in the chair; documents prepared by the general secretary 
(Dr, Ruppel, W, Germany) were first dealt with, 


The most important of these dealt with amendments to the statutes considered at Dusseldorf, introduced by 
by the United Kingdom, U,S,A, and West Germany, That of the first was in the nature of a correction, whereas 
the two others were of a more radical character, The UK delegate proposed to extend the functions of the 
federation to computers, and to include “computation” in the title of the federation, 


The Soviet delegate opposed this proposal, on the grounds that an independent computer federation already 
existed, and that workers in automatic control were only interested in computers insofar as these were component 
parts of control systems, All present agreed with this view, 


The second proposal from the UK delegate amounted to including in the objects of the federation the 
study and development of technical automatic control apparatus, The Soviet delegate supported this proposal, 
Amendments by the U.S, delegate were of a very minor character, since the main proposals put forward earlier 
by Professor Oldenburger had already been agreed and endorsed, A very important extra section, dealing with 
the specialized committees of IFAC and with the rules governing their work, was proposed by the West German 
delegate, 


On the other points the committee agreed to recommend to the General Assembly that: 


1) subscriptions should be 125, 250, 500 and 1000 dollars per year, The subscription chosen was to be 
left to the honor of the member country, and would be published in the press; 


2) Switzerland should be the permanent home of IFAC; 
3) Dusseldorf (W, Germany) should be the permanent location of the secretariat, 


On Professor Oldenburger's proposal, amended in the section on membership, the committee agreed to 
recommend to the General Assembly that the Executive Committee of IFAC should be composed of: 


1) H, Chestnut (U,S,A,), president of IFAC; 
2) A. M, Letov (USSR), first vice-president of IFAC; 
3) V. Broide (France), second vice-president of IFAC; 
4) H, Ruppel (W, Germany), secretary; 
5) Lemann (France), treasurer; 
6) P, Nowatski (Poland); 
1) Gericke (Switzerland); 
8) J, Coals (UK); 
9) Tsien-Hsue-Sen (China); 
10) Evangelisti (Italy); 
11) M, Ainbinder (Belgium), 


This list was later accepted at the General Assembly. 
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Finally, the committee decided to call the first International Congress of IFAC for 1960 in Moscow, 


General Assembly of September 10, 11, and 12, 1957. 


The formal opening of the meeting took place on September 10 in the evening, Representatives of the 
French scientific community welcomed the delegates, The working sessions opened on the following day (the 
11th) with Professor Broide in the chair, A representative of the preparatory committee introduced all the 
delegates and observers, The assembly then considered the foundation statutes proposed by the committee, 
and the draft resolutions drawn up at its earlier meetings, All these resolutions were adopted, 


Then the representative of the General Assembly, Professor Broide, formally declared IFAC in existence, 
and that the statutes of IFAC were submitted for ratification by all participating states, and by all others which 
had either sent observers, or which had indicated their desire to join IFAC prior to September 9, These states 
were Czechoslovakia, Yugoslavia, Chili, Finland, the German Democratic Republic and Denmark, 


The First Meeting of the Executive Committee 


The first meeting of the Executive Committee was held after the close of the General Assembly, 

The following topics were discussed, 

1, The structure of the first IFAC Congress, 

2, Creation of advisory committees on scientific problems by the executive committee, 

3, Time and place of the next meeting of the executive committee, 

Professor Letov presented a preliminary view of the first topic, This was briefly as follows, 

The Congress must cover three aspects of the scientific work of IFAC: 

a) theory and methods of automatic control and regulation; 

b) technical devices of automatic control and related problems; 

c) industrial uses of automatic control and regulation (including the use of computer devices), 

This view was subsequently expanded as to precise scientific content, 

The first resolution on the structure of the Congress must be taken at its next session, set for March, 1958, 
The International Federation of Automatic Control was thus created by resolution of its General Assembly, 


The Soviet scientific community will undoubtedly welcome the creation of IFAC, and is now presented 
with an event of major scientific importance: the first international congress on automatic control in Moscow 
in 1960, 


Specialists in all branches of automatics must begin now to prepare for this congress, in order to demonstrate 
the advanced level of theory and practice in automatic control and regulation in the USSR, and to prepare for 
showing to foreign scientists and engineers everything that is good and useful in our work, 


A, M, Letov and B, N, Naumov. 
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CONFERENCE ON AUTOMATIC CONTROL AND COMPUTER TECHNIQUE 


March 5-8, 1957, a Conference on Automatic Control and Computer Technique called by the All-Union 
Scientific Instrument Engineering Society was held in Moscow, 


Nine hundred delegates from sections of the Academy of Sciences, and from the Academies of the Union 
Republics, higher teaching institutions and research institutes, design offices and laboratories of the Ministries 
and Councils participated, 


Some 40 papers were read at the plenary and sectional meetings, 


The introductory address was given by Rakovskii, who dealt with the general situation in automatic con- 
trol and computer technique in the USSR and abroad and considered detailed examples of the work going on, 


Much attention was devoted to problems related to cybernetics, 


The first part of the address from Solodovnikov (Central Unified Automation Research Institute, as re- 
presentative of the organizing committee) dealt with the tasks and program of the conference; the second dealt 
with defining the object and typical scientific content of technical cybernetics, He observed that the section 
of cybernetics concerned with theory and design principles of machines which effect algoritams can be termed 
technical cybernetics, and that technical cybernetics is composed of the following three fields: 


1) information theory, 


2) theory of logical and computing devices involving studies on methods of transforming information in 
accordance with set algorithms, 


3) theory of automatic control, involving studies on methods of using information for effecting set 
purposes, 


Liapunov demonstrated the importance of cybernetics as being the basis of unified automation, He 
stated that the theory of automatic control played the main part in cybernetics; this in its turn was determined 
by the nature and amount of the circulating information, The flux of circulating information increases with 
the complexity of the control system; optimal information coding problems arise, 


This problem can be dealt with by using the coding methods employed in automatic translation; in parti- 
cular, information should be stored in one code and handled in another, He also pointed out that multiple-step 
automatic control systems have advantages, and gave some examples indicating a need for studies on the control 
processes in living organisms, 


The theory and design principles of control computers were dealt with in papers by Kazakevich, Batkov, 
Kuzin, Shreider, and by Solodovnikov, Batkov, Bredis, and Matveev in a joint paper, 


Kazakevich dealt with "The principles and systems of extremum control,” and presented the main con- 
cepts in the theory of extremum control, He remarked that these systems are fundamental to self-adjusting 
control systems, and gave several examples of peak-holding regulators, 


The paper by Solodonikov, Batkov, Bredis, and Matveev (Central Unified Automation Research Institute) 
dealt with "The current state of the theory of optimal dynamical systems subject to random perturbations;" 
the importance of information theory to control system calculation and design was demonstrated, and the main 
ideas on which statistical methods of studying automatic control systems are based were reviewed, 
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Kuzin dealt with the Z-transform apparatus as used in digital computer automatic control system analysis 
and synthesis, He showed that these operations could thus be performed by methods analogous to those for analog 
computers, An example was considered of a thermal aircraft detector and follower system to illustrate the 
theoretical concepts, 


Batkov dealt with a new application of electronic analogs for determining the correlation function and 
dispersion of the nonstandard random output from an automatic control system from the given characteristics 
of the input perturbation, 


Shreider dealt with the design of "learning" machines, the main feature of which is that they can find 
the optimum control modes from accumulated experience, A necessary feature of such a system is that its 
bela. ior is to a certain extent undefined, 


Dikushin, Liubimirskii, Kamynin and Shtarkman dealt with the development of the most logical automatic 
control systems for production processes, particularly those of metal-working machines, 


Academician Dikushin (Metal-Cutting Machines Experimental Research Institute) dealt with the design of 
program controls, and indicated that proper choice of displacement and speed transducers for the effectors was 
essential, Displacement transducers must require the minimum memory volume as well as being as accurate 
as necessary, The proper form for the speed curves should be one of two steps; in the first the speed is large at 
large displacements, in the second it is less (as the set point is approached), 


Dikushin's idea of "reflex control," the essence of which is that the effect is not controlled from the dis- 
placement direct, but from statistical observations, attracted much attention, 


Liubimirskii, Kamynin, and Shtarkman's report (from the Steklov Mathematics Institute, Academy of 
Sciences of the USSR) dealt with the optimal coding of information in automatic control systems (including 
multiposition process controls), It is shown that for a lathe, for instance, doing some simple operation, no 
complex control was needed; the functions could be performed by a simple device, However, to supervise the 
functioning of the device a rather complex universal control machine is needed, but this can serve for several 
objects, and so the system is much simplified, 


In conclusion Liubimirskii observed that the algorithms now available (e.g, for machine translation, 
for solving mathematical problems, etc,) are much more complex that those needed for many industrial pro- 
cesses; he pointed out that joint teams of mathematicians, technologists and computer specialists should be set 
up to solve process algorithm problems rapidly, and to design control computers, 


The uses of high-speed computers are not limited to automatic process control or to mathematical problems, 
Computers can be used to process any information in any way, including translations from one language to another, 
programming programs, etc, and extensive development prospects are to be envisaged, 


Shur-Bur's paper (from the Steklov Mathematics Institute, Academy of Sciences of the USSR) dealt with 
these topics, 


A second extensive group of papers dealt mainly with the components and common features of computers, 
with programming methods, and with calculations and problems encountered in the automatic control field, 


Recent advances in producing high-speed, reliable, cheap small computer components were dealt with by 
Korol'kov, by Mamonovy and by Sharapov, 


The first deal with some very promising high-speed magnetic switching elements using rectangular hysteresis 
loop materials, It was also pointed out that materials not having rectangular hysteresis loops could be used, 
Standard block units performing arithmetical, logical, etc, operations could be used in automatic devices, 


The second dealt with semiconductor devices used in digital computers, The special features of point-contact 
and planar germanium transistors operating at large signals, and the various designs of trigger, were considered, 
Transistor triggers using carrier storage in planar diodes and the S-shaped response curves of point-contact diodes, 
circuits of units comprised of point-contact transistor dynamic triggers, or planar-contact static triggers (recorders, 
shift registers, read-out units) were presented, 


Zimin's and Guienmakher's papers dealt with the use of these components in computer circuits, 
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Gutenmakher showed that magnetic switching units can provide long-term and intermediate storage 
facilities, arithmetical, logical, etc, operations, or control and switching operations, Various circuits were 
illustrated for adding n-column binary numbers, for multipliers of series- parallel action for two n- column 
binary numbers, circulating registers, decimal recorders, switching devices, logical units, etc, all using 
magnetic switching units, 


Zimin dealt with semiconductor devices used in computers, and pointed out that scientifically sound 
specifications for such devices should be developed, 


Dealing with machine logic, he showed that all logical operations the machine can perform can be built 
up from the three basic elements “or,” "and," and "no," These elements can be made to give low attenuation 
and satisfactory signal-to-noise ratio using diode plus transformer circuits, Registers and recorders using transistor 
amplifiers can be built up on these principles, 


The input and output units of computers were considered by Ryzhov, by Trubnikov and by Zavolokin and 
Baskukov, 


Ryzhov reviewed the foreign published data on analog-digital and digital-analog converters, 


Trubnikov dealt with the comparative aspects of the input and output units of high-speed program-con- 
trolled computers, and analyzed the various forms of program carrier, 


Zavolokin and Baskukov dealt with the input and output units of cyclic converters with inductive phase- 
shifters, 


Val'denberg dealt with a specialized analog computer specially designed to solve the Fredholm and 
Volterra integral equations of first and second orders with nuclei of configuration type so often encountered 
in control problems, 


The iterative Gauss-Seidel method was used, Digital elements were widely used, The machine found 
the solution periodically, without preserving the natural time-scale, 


Mathematical statistics techniques can be used tc ° ~' with many physical phenomena, However, the 
great labor involved in computing many statistical parat. _ts of random processes (particularly correlation 
functions) prevents their wide use in practice, 


Novikov (Institute of Automation and Remote Control, Academy of Sciences of the USSR) has produced 
a specialized computer (a magnetic correlograph) for computing correlation functions automatically; the 
design of the device was described, The signal was recorded on a magnetic tape, covering the frequency 
range from dc up to 200 cps. The cycling time varied from 10 to 60 minutes, 


Vitenberg described an automatic electronic analog for seeking the solution to a given system of 
equations, In relation to computers used as optimum solution seekers, he stated that electronic analogs were 
the best for use in such cases, 


Survey and minimalization methods are the best for finding optimum solutions, An electronic analog 
for seeking the solution to a given system of equations consists of the analog proper plus control system and 
set of initial conditions, a timing unit, and a unit for controlling the logical operations of the input and output 
devices, 


The comparative features and design methods of such devices were considered, and possible ways of pro- 
duction indicated, 


Maiorov and Zhidkov dealt with the mathematical aspects of digital differential analyzers and with 
their use as control devices, 


Zhidkov described how nvmbers were represented in digital differential analyzers and with ways of 
setting up the problems and choosing scales, 


Maiorov dealt with how such analyzers are used in control systems covering many machines, how they 
are coupled to real objects, with their possible uses, and with how they compare with other digital computers, 


Gutenmakher dealt with the prospects of using a new type of information and statistical machine in 
control systems, The rapid development of various types of storé was pointed out, and the structures of future 
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types of information and statistical machine were considered in detail; it was pointed out that new types of 
reading device for such machines will take the lines of analogs to mental processes, 


The applications of future types of information and statistical machine were considered (in chemistry, 
in statistics, in trade) to round off the paper. 


Several papers dealt with computers designed to control actual industrial processes (those of Kopai-Gora, 
Khetagurov, Kaganov, Fitsner and Vasil’ev, Efroimovich and of Cheliustkin), 


Efroimovich (Central Automatics Laboratory), Cheliustkin (Central Automatics Laboratory), and Kopai- 
Gora dealt with computers designed to control the main objects in ferrous metallurgy (blast-furnace processes, 
arc steel- melting furnaces, rolling- mills), 


In all these papers the good prospects for using computers in ferrous metallurgy were pointed out, the 
state of the art abroad was summarized, and the situation in the USSR was shown to be such that the time was 
ripe for introducing computer control in this field, The lack of mathematical descriptions for the processes 
was indicated as greatly retarding this development, 


Several examples of computer devices at present under test were mentioned, Much attention was de- 
voted to future uses of computers, The 1956 All-Union Conference on Process Automation in Ferrous Metallurgy 
(in Magnitogorsk), at which a plan to provide full automatic working for several processes was adopted, was 
shown to have been very important, 


The choice of most suitable cutting conditions is of great importance in increasing productivity in metal- 
working processes, 


Vasil'ev and Fitsner dealt with a computer designed to handle this problem, Up to 20 distinct quantities 
which determine the cutting conditions can be fed in, If some of these are given, the required optimal para- 
meters can be calculated in 2-3 minutes, Laboratory and works test results were stated to be satisfactory, 


Khetaguroy dealt with a digital control system for a metal- working machine, 
After thorough discussions on the papers the Conference adopted several resolutions, 


In view of the primary importance of fully automatic process control, it was resolved to extend research 
on the use of computers in automatic process control, 


In doing this the main attention must be devoted to developing technical cybernetics theory, to producing 
new control machine components and standard units and racks, to developing new process techniques in various 
industries so that computer technique can be used, 


It was also resolved that the Central Unified Automation Research Institute, the research institues, etc, of 
the Ministries and Councils, and the Academy of Sciences, should organize specialist groups on technical cyber- 
netics, It was decided to ask the All-Union Scientific Engineering Society to hold a conference on cybernetics 
in collaboration with the Academy of Sciences in 1958, 


This conference advanced the development of automatic control theory, made available experience with 
particular control problems involving computer techniques, and enabled specialists in different aspects to combine 
forces in order to solve the problems of full-scale process automation in the USSR, 


I, M, Rusevich, 


ERRATA 


In the article by Ia, Z, Tsypkin in Volume 17, No, 12 (1956): 
1, In Formula (2,13) and all following, A“ue[n] should be replaced by (—1)*A*us [n — k] (k = 0, 1,2, -- -)- 


2. Formula (2,19) should read: 
forr=1 

Ujn st =don— dy, 
forr=2 


Us gt = don*— dy (Qn— 1) + 


